CHAPTER 

6 


POWER FLOW ANALYSIS 


6.1 INTRODUCTION 

In the previous chapters, modeling of the major components of an electric power 
system was discussed. This chapter deals with the steady-state analysis of an in- 
terconnected power system during normal operation. The system is assumed to be 
operating under balanced condition and is represented by a single-phase network. 
The network contains hundreds of nodes and branches with impedances specified 
in per unit on a common MVA base. 

Network equations can be formulated systematically in a variety of forms. 
However, the node-voltage method, which is the most suitable form for many 
power system analyses, is commonly used. The formulation of the network equa- 
tions in the nodal admittance form results in complex linear simultaneous algebraic 
equations in terms of node currents. When node currents are specified, the set of 
linear equations can be solved for the node voltages. However, in a power system, 
powers are known rather than currents. Thus, the resulting equations in terms of 
power, known as the power flow equation, become nonlinear and must be solved 
by iterative techniques. Power flow studies, commonly referred to as load flow, are 
the backbone of power system analysis and design. They are necessary for plan- 
ning, operation, economic scheduling and exchange of power between Utilities. In 
addition, power flow analysis is required for many other analyses such as transient 
stability and contingency studies. 
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In this chapter, the bus admittance matrix of the node-voltage equation is 
formulated, and a MATLAB function named ybus is developed for the system- 
atic formation of the bus admittance matrix. Next, two commonly used iterative 
techniques, namely Gauss-Seidel and Newton-Raphson methods for the solution 
of nonlinear algebraic equations, are discussed. These techniques are employed in 
the solution of power flow problems. Three programs Ifgauss, Ifnewton, and de- 
couple are developed for the solution of power flow problems by Gauss-Seidel, 
Newton-Raphson, and the fast decoupled power flow, respectively. 


6.2 BUS ADMITTANCE MATRIX 

In order to obtain the node-voltage equations, consider the simple power system 
shown in Figure 6.1 where impedances are expressed in per unit on a common 
MVA base and for simplicity resistances are neglected. Since the nodal solution is 
based upon Kirchhoff’s current law, impedances are converted to admittance, i.e., 

1 1 

Vij — ~ : : 

Zij Tij + JXij 



FIGURE 6.1 

The impedance diagram of a simple system. 
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FIGURE 6.2 

The admittance diagram for system of Figure 6.1. 


The Circuit has been redrawn in Figure 6.2 in terms of admittances and trans- 
formation to current sources. Node 0 (which is normally ground) is taken as refer- 
ence. Applying KCL to the independent nodes 1 through 4 results in 

h = 2 /ioVi + yn{Vi - V 2 ) + yn{V\ - V 3 ) 
h = 2/20 V2 + 2/i2(Vi — Vi) + 2/23 (V2 — V3) 

0 = 2/23 (V3 — V2) -I- 2 /i 3 (V 3 — Vi) + 2/34 (V3 — V4) 

0 = 2/34 (V4 — V3) 

Rearranging these equations yields 

h = (2/10 + 2/12 + 2/i3)Vi — j/12 V-2 — 2/13V3 
h = — 2/i2Vi + (2/20 + 2/12 + 2/23) V2 — 2 / 23 V3 

0 = —2/13 Vi - 2/23 Vi + (2/13 + 2/23 + 2/34) V3 - 2/34 Vi 
0 = — 2/34V3 + 2/34 V4 

We introduce the following admittances 


Vii = 2/io + 2/12 + 2/13 
V22 = 2/20 + 2/12 + 2/23 
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*33 = 2/13 + 2/23 + 2/34 
*44 = 2/34 
*12 = *21 = “ 2/12 
*"13 = *31 = —2/13 
*23 = *32 = —2/23 
*34 = *43 = 2/34 

The node equation reduces to 

h = *11*1 + *12*2 + *13*3 + *14*4 
I 2 = *21*1 + *22*2 + *23*3 + *24*4 
h = *31*1 + *32*2 + *33*3 + *34*4 
I 4 — *41*1 + *42*2 + *43*3 + *44*4 

In the above network, since there is no connection between bus 1 and 4, *1 4 = 
*41 = 0; similarly *24 = *42 = 0. 

Extending the above relation to an n bus system, the node-voltage equation 
in matrix form is 

' h 
I 2 

li 

. J n 
or 


where Ifc us is the vector of the injected bus currents (i.e., extemal current sources). 
The current is positive when flowing towards the bus, and it is negative if flowing 
away from the bus. V^ us is the vector of bus voltages measured from the reference 
node (i.e., node voltages). Yf, us is known as the bus admittance matrix. The diag- 
onal element of each node is the sum of admittances connected to it. It is known as 
the self-admittance or driving point admittance, i.e., 

n 

*« = Vij i ^ * (6.3) 

3=0 

The off-diagonal element is equal to the negative of the admittance between the 
nodes. It is known as the mutual admittance or transfer admittance, i.e., 
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( 6 . 1 ) 


lims — V fnis 


( 6 . 2 ) 


*1? “ *ji — Vij 


( 6 . 4 ) 
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When the bus currents are known, (6.2) can be solved for the n bus voltages. 

V blis = Ybus I bus (6.5) 


The inverse of the bus admittance matrix is known as the bus impedance matrix 
Zbus- The admittance matrix obtained with one of the buses as reference is nonsin- 
gular. Otherwise the nodal matrix is singular. 

Inspection of the bus admittance matrix reveals that the matrix is symmetric 
along the leading diagonal, and we need to store the upper triangular nodal ad- 
mittance matrix only. In a typical power system network, each bus is connected to 
only a few nearby buses. Consequently, many off-diagonal elements are zero. Such 
a matrix is called sparse, and efficient numerical techniques can be applied to com- 
pute its inverse. By means of an appropriately ordered triangular decomposition, 
the inverse of a sparse matrix can be expressed as a product of sparse matrix fac- 
tors, thereby giving an advantage in computational speed, storage and reduction of 
round-off errors. However, Z^us-, which is required for short-circuit analysis, can be 
obtained directly by the method of building algorithm without the need for matrix 
inversion. This technique is discussed in Chapter 9. 

Based on (6.3) and (6.4), the bus admittance matrix for the network in Figure 
6.2 obtained by inspection is 


Ybus 


-j 8.50 y'2.50 j5.00 

j2.50 — j'8.75 j5.00 

j5.00 j5.00 — j'22.50 

0 0 y‘12.50 


0 ' 
0 

y'12.50 
— j'12.50 . 


A function called Y = ybus(zdata) is written for the formation of the bus 
admittance matrix, zdata is the line data input and contains four columns. The 
first two columns are the line bus numbers and the remaining columns contain the 
line resistance and reactance in per unit. The function retums the bus admittance 
matrix. The algorithm for the bus admittance program is very simple and basic to 
power system programming. Therefore, it is presented here for the reader to study 
and understand the method of solution. In the program, the line impedances are 
first converted to admittances. Y is then initialized to zero. In the first loop, the 
line data is searched, and the off-diagonal elements are entered. Finally, in a nested 
loop, line data is searched to find the elements connected to a bus, and the diagonal 
elements are thus formed. 

The following is a program for building the bus admittance matrix: 
function [Y] = ybus (zdata) 

nl=zdata( : , 1) ; nr=zdata( : ,2) ; R=zdata( : ,3) ; X=zdata( : ,4) ; 
nbr=length (zdata ( : ,1)) ; nbus = max (max (ni) , max(nr)); 

Z = R + j*X; “/.branch impedance 
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y= onesfnbr ,1) ./Z; 7,branch admittance 

Y = zeros(nbus,nbus) ; 7, initialize Y to zero 

for k = l:nbr; 7« formation of the off diagonal elements 
if nl(k) > 0 fe nr(k) > 0 
Y(nl(k) ,nr(k)) = Y(nl(k) ,nr(k)) - y(k); 

Y(nr(k) ,nl(k)) = Y(nl(k) ,nr(k) ) ; 
end 

end 

for n = l:nbus 7. formation of the diagonal elements 

for k = 1 :nbr 
if nl(k) == n | nr(k) == n 
Y(n,n) = Y(n,n) + y(k) ; 
else, end 
end 

end 

Example 6.1 

The emfs shown in Figure 6.1 are E\ = 1.1/0° and — 1.0Z0 0 . Use the func- 
tion Y = ybus(zdata) to obtain the bus admittance matrix. Find the bus impedance 
matrix by inversion, and solve for the bus voltages. 


With source transformation, the equivalent current sources are 


h 

h-- 


1.1 
= yl.O 
1.0 

708 : 


= -j 1.1 pu 
— y'1.25 pu 


The following commands 

7. From To R X 
z = [ 0 1 0 1.0 

0 2 0 0.8 

1 2 0 0.4 

1 3 0 0.2 

2 3 0 0.2 

3 4 0 0.08]; 

Y = ybus(z) 7. bus admittance matrix 

Ibus = C— j *1.1; -j*1.25; 0; 0]; 7« vector of bus currents 

Zbus = inv(Y) 7. bus impedance matrix 

Vbus = Zbus* Ibus 


result in 
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Y = 


0 

- 

8 . 50i 

0 

+ 

2. 50i 

0 

+ 

5. OOi 

0 

+ 

O.OOi 

0 

+ 

2 . 50i 

0 

- 

8 . 75i 

0 

+ 

5. OOi 

0 

+ 

O.OOi 

0 

+ 

5 . OOi 

0 

+ 

5. OOi 

0 

- 

22 . 50i 

0 

+ 

12. 50i 

0 

+ 

O.OOi 

0 

+ 

O.OOi 

0 

+ 

12. 50i 

0 

- 

12. 50i 

Zbus = 












0 

+ 

0 . 50i 

0 

+ 

0.40i 

0 

+ 

0 . 450i 

0 

+ 

0.450i 

0 

+ 

0.40i 

0 

+ 

0 . 48i 

0 

+ 

0 . 440i 

0 

+ 

0 . 440i 

0 

+ 

0.45i 

0 

+ 

0 . 44i 

0 

+ 

0 . 545i 

0 

+ 

0 . 545i 

0 

+ 

0.45i 

0 

+ 

0 . 44i 

0 

+ 

0 . 545i 

0 

+ 

0 . 625i 


Vbus = 

1.0500 
1 . 0400 
1 . 0450 
1 . 0450 

The solution of equation I^s = ^bus^bus by inversion is very inefficient. It 
is not necessary to obtain the inverse of Y^. Instead, direct solution is obtained 
by optimally ordered triangular factorization. In MATLAB, the solution of linear 
simultaneous equations AX = B is obtained by using the matrix division operator 
\ (i.e., X = A \ B), which is based on the triangular factorization and Gaussian 
elimination. This technique is superior in both execution time and numerical accu- 
racy. It is two to three times as fast and produces residuals on the order of machine 
accuracy. 


In Example 6.1, obtain the direct solution by replacing the statements Zbus = 
inv(Y) and Vbus = Zbus*Ibus with Vbus = Y\ Ibus. 


6.3 SOLUTION OF NONLINEAR 
ALGEBRAIC EQUATIONS 

The most common techniques used for the iterative solution of nonlinear algebraic 
equations are Gauss-Seidel, Newton-Raphson, and Quasi-Newton methods. The 
Gauss-Seidel and Newton-Raphson methods are discussed for one-dimensional 
equation, and are then extended to n-dimensional equations. 

6.3.1 GAUSS-SEIDEL METHOD 

The Gauss-Seidel method is also known as the method of successive displace- 
ments. To illustrate the technique, consider the solution of the nonlinear equation 
given by 


f(x) = 0 


( 6 . 6 ) 
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The above function is rearranged and written as 

x = g(x ) (6.7) 

If is an initial estimate of the variable x, the following iterative sequence is 
formed. 


x (fc+1 ) = g(x {k) ) (6.8) 

A solution is obtained when the difference between the absolute value of the suc- 
cessive iteration is less than a specified accuracy, i.e., 

|se ( * +1 > - *<*>| < c (6.9) 

where e is the desired accuracy. 


Example 6.2 

Use the Gauss-Seidel method to find a root of the following equation 
f(x) = x 3 — 6x 2 + 9x — 4 = 0 


Solving for x, the above expression is written as 


1 3 6 2 4 

* = “9* + 9* + 9 


= «(*) 


The MATLAB plot command is used to plot g(x) and x over a range of 0 to 4.5, 
as shown in Figure 6.3. The intersections of g(x) and x results in the roots of 
/(x). From Figure 6.3 two of the roots are found to be 1 and 4. Actually, there 
is a repeated root at x = 1. Apply the Gauss-Seidel algorithm, and use an initial 
estimate of 


x<°> = 2 

From (6.8), the first iteration is 

x^ = 5 (2) = -i(2) 3 + jj(2) 2 + ^= 2.2222 

The second iteration is 

x (2 > = 3(2.2222) = -J(2.2222) 3 + ^(2.2222) 2 + £ = 2.5173 

y y y 

The subsequent iterations result in 2.8966, 3.3376, 3.7398, 3.9568, 3.9988 and 
4.0000. The process is repeated until the change in variable is within the desired 
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X 


FIGURE 6.3 

Graphical illustration of the Gauss-Seidel method. 


accuracy. It can be seen that the Gauss-Seidel method needs many iterations to 
achieve the desired accuracy, and there is no guarantee for the convergence. In this 
example, since the initial estimate was within a “boxed in” region, the solution 
converged in a zigzag fashion to one of the roots. In faet, if the initial estimate 
was outside this region, say x^ = 6, the process would diverge. A test of conver- 
gence, especially for the n-dimensional case, is difficult, and no general methods 
are known. 

The following commands show the procedure for the solution of the given 
equation starting with an initial estimate of x^ = 2 . 


7. Change in variable is set to a high value 

’/, Initial estimate 
7. Iteration counter 
x’)7iHeading for results 
< 100 7, Test for convergence 
7. No. of iterations 


g 

>= 0.001 
i; 


dx 

iter 


dx=l ; 
x=2; 

iter = 0; 
disp( ’ Iter 
while abs(dx) 
iter = iter + 
g = -l/9*x~3+6/9*x~2+4/9 ; 
dx = g-x; 
x = x + dx; 

f printf ( ’7.g’ , iter), disp([g, 
end 


7. Change in variable 
7« Successive approximation 
dx, x]) 


The result is 
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Iter 

g 

dx 

X 

1 

2.2222 

0 . 2222 

2.2222 

2 

2.5173 

0.2951 

2.5173 

3 

2.8966 

0 . 3793 

2.8966 

4 

3 . 3376 

0.4410 

3.3376 

5 

3.7398 

0.4022 

3.7398 

6 

3.9568 

0.2170 

3.9568 

7 

3.9988 

0 . 0420 

3.9988 

8 

4.0000 

0.0012 

4.0000 

9 

4.0000 

0.0000 

4.0000 


In some cases, an acceleration factor can be used to improve the rate of conver- 
gence. If a > 1 is the acceleration factor, the Gauss-Seidel algorithm becomes 


x (fc+l) _ x (k) a [g( X ( k )} _ a;( fc )] (6.10) 

Example 6.3 

Find a root of the equation in Example 6.2, using the Gauss-Seidel method with an 
acceleration factor of a = 1.25: 


Starting with an initial estimate of = 2 and using (6.10), the first iteration is 

g ( 2 ) = -|( 2) 3 + |( 2) 2 + | = 2.2222 

x (1) = 2 + 1.25[2.2222 - 2] = 2.2778 

The second iteration is 

16 4 

g( 2.2778) = --(2.2778) 3 + -(2.2778) 2 + - = 2.5902 

x (2) = 2.2778 + 1.25[2.5902 - 2.2778] = 2.6683 


The subsequent iterations result in 3.0801, 3.1831, 3.7238, 4.0084, 3.9978 and 
4.0005. The effect of acceleration is shown graphically in Figure 6.4. Care must 
be taken not to use a very large acceleration factor since the larger step size may 
result in an overshoot. This can cause an increase in the number of iterations or 
even result in divergence. In the MATLAB command of Example 6.2, replace the 
command before the end statement by x = x + 1.25 * dx to reflect the effect of the 
acceleration factor and run the program. 
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FIGURE 6.4 

Graphical illustration of the Gauss-Seidel method using acceleration factor. 

We now consider the system of n equations in n variables 

,x n ) = Cl 

/ 2 (xi,x 2 , •••,£„) = c 2 (6.11) 


fn(x 1 j X 2) ’ ‘ ' , X n ) — C n 

Solving for one variable from each equation, the above functions are rearranged 
and written as 

Xl = Cl + ffl (xi, X 2 , ' ' ' j X n ) 

^2 = C 2 + g2(xi,X2, • • ■ , x n ) (6.12) 


— Cn ~ H 5n(^l) 3-2> ' ' ' i Xn) 

The iteration procedure is initiated by assuming an approximate solution for each 
of the independent variables (x{°\ • • ■ ,x^). Equation (6.12) results in a new 

approximate solution (xf \ xjj 1 ^ • • • , x^). In the Gauss-Seidel method, the updated 
values of the variables calculated in the preceding equations are immediately used 
in the solution of the subsequent equations. At the end of each iteration, the cal- 
culated values of all variables are tested against the previous values. If all changes 





200 6. POWER FLOW ANALYSIS 


in the variables are within the specified accuracy, a solution has converged, oth- 
erwise another iteration must be performed. The rate of convergence can often be 
increased by using a suitable acceleration factor a, and the iterative sequence be- 
comes 




= x- 




(6.13) 


6.3.2 NEWTON-RAPHSON METHOD 

The most widely used method for solving simultaneous nonlinear algebraic equa- 
tions is the Newton-Raphson method. Newton’s method is a successive approxima- 
tion procedure based on an initial estimate of the unknown and the use of Taylor’s 
series expansion. Consider the solution of the one-dimensional equation given by 

f{x ) = c (6.14) 

If x^ 0 ' is an initial estimate of the solution, and Ax^ is a small deviation from the 
correct solution, we must have 

f{x^ + Ax (°)) = c 

Expanding the left-hand side of the above equation in Taylor’s series about x- 0 ^ 
yields 

i (g) V-)-...... 

Assuming the error Ax^ 0 ) is very small, the higher-order terms can be neglected, 
which results in 



where 

Ac<°> =c-/(x ( °)) 

Adding Ax^ to the initial estimate will result in the second approximation 

A c (°) 

(dL\ i0) 


= x (°) + 
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Successive use of this procedure yields the Newton-Raphson algorithm 

Ac (fc) = c - /(x (fc) ) (6.15) 


Ax< fe) = 


AcW 



(6.16) 


x(fc+!) _ x (k) + Ax (k) 

(6.17) 

(6.16) can be rearranged as 


Ac ( fc ) = j( fc )Ax (fc) 

(6.18) 

where 




The relation in (6.18) demonstrates that the nonlinear equation f(x) — c = 0 is 
approximated by the tangent line on the curve at x^. Therefore, a linear equation 
is obtained in terms of the small changes in the variable. The intersection of the 
tangent line with the x-axis results in x^ fc+1 ). This idea is demonstrated graphically 
in Example 6.4. 

Example 6.4 

Use the Newton-Raphson method to find a root of the equation given in Example 
6.2. Assume an initial estimate of x^ = 6. 


The MATLAB plot command is used to plot f(x ) = x 3 — 6x 2 + 9x — 4 over 
a range of 0 to 6 as shown in Figure 6.5. The intersections of f(x) with the x-axis 
results in the roots of /(x). From Figure 6.5, two of the roots are found to be 1 and 
4. Actually, there is a repeated root at x = 1. 

Also, Figure 6.5 gives a graphical description of the Newton-Raphson method. 
Starting with an initial estimate of x^ = 6, we extrapolate along the tangent to 
its intersection with the x-axis and take that as the next approximation. This is 
continued until successive x-values are sufficiently close. 

The analytical solution given by the Newton-Raphson algorithm is 


df (g) 

dx 


= 3x 2 


12x + 9 


Ac (0) = c - /(x (0) ) = 0 - f(6) 3 - 6(6) 2 + 9(6) - 4] = -50 
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FIGURE 6.5 

Graphical illustration of the Newton-Raphson algorithm. 



AxW 


Ac<°> _ -50 

(4O (0) 


45 


= - 1.1111 


J x ) 


Therefore, the result at the end of the first iteration is 

x {1) = z<°) + Ax (0) = 6 — 1.1111 = 4.8889 
The subsequent iterations result in 

æ (2) _ x (0 + Ax^ = 4.8889 - = 4.2789 

22.037 


x (3) = x (2) + Ax (2) = 4.2789 - = 4.0405 

12.5797 


Æ ( 4 ) _ x (3) _j_ Ax^ = 4.0405 — ^ = 4.0011 

9.4914 

x (5) = x (4) + Ax (4) = 4.0011 - = 4.0000 

9.0126 
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We see that Newton’s method converges considerably more rapidly than the Gauss- 
Seidel method. The method may converge to a root different from the expected one 
or diverge if the starting value is not close enough to the root. 


The following commands show the procedure for the solution of the given equation 
by the Newton-Raphson method. 

dx=l; 7, Change in variable is set to a high value 

x=input( ’Enter initial estimate -> ’); 7. Initial estimate 
iter = 0; 7. Iteration counter 

disp(’iter Dc J dx x’) 7. Heading 

while abs(dx) >= 0.001 & iter < 100 Test for convergence 
iter = iter +1; 7. No. of iterations 

Dc = 0 - (x“3 - 6*x~2 + 9*x - 4) ; '/, Residual 

J = 3*x~2 - 12*x + 9; 7. Derivative 

dx= Dc/J; 7.Change in variable 

x=x + dx; 7. Successive solution 

f printf ( ’7.g’ , iter), disp([Dc, J, dx, x] ) 
end 

The result is 


Enter the initial estimate -> 6 


iter 

Dc 

J 

dx 

X 

1 ■ 

-50.0000 

45 . 0000 

- 1.1111 

4.8889 

2 ■ 

-13.4431 

22 . 0370 

-0.6100 

4.2789 

3 

-2.9981 

12.5797 

-0.2383 

4.0405 

4 

-0.3748 

9.4914 

-0.0395 

4.0011 

5 

-0.0095 

9.0126 

-0.0011 

4.0000 

6 

- 0.0000 

9.0000 

- 0.0000 

4 . 0000 


Now consider the n-dimensional equations given by (6.11). Expanding the left- 
hand side of the equations (6.11) in the Taylor’s series about the initial estimates 
and neglecting all higher order terms, leads to the expression 

(/!)( o, + (g)%,<o, + (g)%4o. + ... + (g)% 4 , = C2 


204 6. POWER FLOW ANALYSIS 


or in matrix form 


’ Cl - (/l) (0) ' 


’ (åh} 

(0) 

.. 

. (åh} 

(0)1 





\9a:i j 



\dXnj 




C2 - (/2) (0) 

— 

(åh} 

r 

1 9X2 J 

' 

l (0) 


a4 0) 

. C„ ~ (/n) (0) . 


_ (åk, 

l (0) 

(fef " 


(0) 




In short form, it can be written as 

A C (fc) = J (fc) AX (fe) 


or 

AX (fc) = [jW] _1 AC7 (fc) 


(6.19) 


and the Newton-Raphson algorithm for the n-dimensional case becomes 
where 



Ax ^ 


ci - (/i) (fc) 

AX^ = 

Ax? 

and A C {k) = 

C2 - (/2) (fc) 


AXn ' 1 


. Cn - (/„)<*> . 


( 6 . 20 ) 


( 6 . 21 ) 



( 6 . 22 ) 


jr(fc) 

is called the Jacobian matrix. Elements of this matrix are the partial 
derivatives evaluated at X^ k \ It is assumed that j( fc ) has an inverse during each 
iteration. Newton’s method, as applied to a set of nonlinear equations, reduces the 
problem to solving a set of linear equations in order to determine the values that 
improve the accuracy of the estimates. 

The solution of (6.19) by inversion is very inefficient. It is not necessary 
to obtain the inverse of j( k \ Instead, a direct solution is obtained by optimally 
ordered triangular factorization. In MATLAB, the solution of linear simultaneous 
equations A C = JAX is obtained by using the matrix division operator \ (i.e., 
AX = .J \ AG') which is based on the triangular factorization and Gaussian elim- 
ination. 
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Example 6.5 

Use the Newton-Raphson method to find the intersections of the curves 

x\ + x 2 2 = 4 
e Xl + x 2 = l 

Graphically, the solution to this system is represented by the intersections of the 
circle x\ + x\ = 4 with the curve e Xl + x 2 — 1. Figure 6.6 shows that these are 
near (1, —1.7) and (—1.8, 0.8). 



FIGURE 6.6 
Graphs of Example 6.5. 


Taking partial derivatives of the above functions results in the Jacobian matrix 


J 


2x\ 2 x 2 


o x l 


The Newton-Raphson algorithm for the above system is presented in the following 
statements. 


iter = 0; 


*/. Iteration counter 
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x=input( , Enter initial estimates, col. vector [xl ;x2] ->’ ) ; 
Dx = [1; 1] ; 7« Change in variable is set to a high value 

C= [4; 1]; 

dispC’Iter DC Jacobian matrix Dx x’); 

7, Heading for results 

while max(abs(Dx)) >= 0.0001 & iter <10 7,Convergence test 
iter=iter+l; 7, Iteration counter 

f = [x(l)~2+x(2)~2; exp(x(l) )+x(2)] ; 7. Functions 

DC = C - f ; 7. Residuals 

J = [2*x(l) 2*x(2) 7. Jacobian matrix 

exp(x(l)) 1]; 

Dx=J\DC; 7. Change in variables 

x=x+Dx ; 7. Successive solutions 

f printf ( ’7*g’ , iter), disp([DC, J, Dx, x]) 7. Results 

end 


When the program is run, the user is prompted to enter the initial estimate. 
Let us try an initial estimate given by [0.5; -1], 


Enter Initial estimates, col. vector [x \ ; X 2 ] — * [0.5; — 1] 


Iter 

A C 

Jacobian 

matrix 

Ax 

X 

1 

2.7500 

1 . 0000 

-2.0000 

0 . 8034 

1 . 3034 


0.3513 

1 . 6487 

1.0000 

-0.9733 

-1.9733 

2 

-1.5928 

2 . 6068 

-3.9466 

-0.2561 

1 . 0473 


-0.7085 

3.6818 

1 . 0000 

0 . 2344 

-1.7389 

3 

-0.1205 

2.0946 

-3.4778 

-0.0422 

1.0051 


- 0.1111 

2.8499 

1.0000 

0.0092 

-1.7296 

4 

-0.0019 

2.0102 

-3.4593 

-0.0009 

1 . 0042 


-0.0025 

2.7321 

1 . 0000 

0.0000 

-1.7296 

5 

- 0.0000 

2.0083 

-3.4593 

- 0.0000 

1 . 0042 


- 0.0000 

2.7296 

1.0000 

- 0.0000 

-1 . 7296 


After five iterations, the solution converges to x\ = 1.0042 and = —1.7296 
accurate to four decimal places. Starting with an initial value of [—0.5; 1], which 
is doser to the other intersection, results in x\ = —1.8163 and x% — 0.8374. 

Example 6.6 

Starting with the initial values, x\ = 1, X 2 = 1, and x 3 = 1, solve the following 
system of equations by the Newton-Raphson method. 

x\ — x\ + xl = 11 
X\X 2 + x\ — 3^3 = 3 

xi - X1X3 + X2X3 = 6 
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Taking partial derivatives of the above functions results in the Jacobian matrix 


J = 


2x\ — 2 x 2 2 x 3 

x 2 xi + 2 x 2 —3 

1 - X 3 x 3 X\ + x 2 


The following statements solve the given system of equations by the Newton- 
Raphson algorithm 


Dx=[10; 10; 10] ; y.Change in variable is set to a high value 
x=[l; 1; 1]; % Initial estimate 

C« [11; 3; 6]; 

iter = 0; '/, Iteration counter 

while max (abs(Dx))>=. 0001 & iter<10 ;’/,Test for convergence 
iter = iter + 1 ’/, No. of iterations 

F = [x(l)~2-x(2)“2+x(3)''2 */, Functions 

x(l)*x(2)+x(2) ~2-3*x(3) 
x(l)-x(l)*x(3)+x(2)*x(3)] ; 


DC =C - F 
J = [2*x(l) 
x(2) 
l-x(3) 
Dx=J\DC 
x=x+Dx 


7. Res i dual s 

-2*x(2) 2*x(3) 7, Jacobian matrix 

x(l)+2*x(2) -3 

x(3) -x(l)+x(2)] 

7,Change in variable 
7. Successive solution 


end 


The program results for the first iteration are 


DC = 

10 

4 

5 

Dx = 

4.750 

5.000 

5.250 


2-2 2 
13-3 
0 10 

5.750 

6.000 

6.250 


After six iterations, the solution converges to x\ — 2.0000, x 2 = 3.0000, and 
x 3 = 4.0000. 

Newton’s method has the advantage of converging quadratically when we 
are near a root. However, more functional evaluations are required during each 
iteration. A very important limitation is that it does not generally converge to a 
solution from an arbitrary starting point. 
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6.4 POWER FLOW SOLUTION 

Power flow studies, commonly known as load flow, form an important part of 
power system analysis. They are necessary for planning, economic scheduling, and 
control of an existing system as well as planning its future expansion. The problem 
consists of determining the magnitudes and phase angle of voltages at each bus and 
active and reactive power flow in each line. 

In solving a power flow problem, the system is assumed to be operating under 
balanced conditions and a single-phase model is used. Four quantities are associ- 
ated with each bus. These are voltage magnitude |F|, phase angle S, real power P, 
and reactive power Q. The system buses are generally classified into three types. 


Slack bus One bus, known as slack or swing bus, is taken as reference where the 
magnitude and phase angle of the voltage are specified. This bus makes up 
the difference between the scheduled loads and generated power that are 
caused by the losses in the network. 

Load buses At these buses the active and reactive powers are specified. The mag- 
nitude and the phase angle of the bus voltages are unknown. These buses are 
called P-Q buses. 

Regulated buses These buses are the generator buses. They are also known as 
voltage-controlled buses. At these buses, the real power and voltage magni- 
tude are specified. The phase angles of the voltages and the reactive power 
are to be determined. The limits on the value of the reactive power are also 
specified. These buses are called P-V buses. 


6.4.1 POWER FLOW EQUATION 

Consider a typical bus of a power system network as shown in Figure 6.7. Trans- 
mission lines are represented by their equivalent 7r models where impedances have 
been converted to per unit admittances on a common MVA base. 

Application of KCL to this bus results in 


h — ViC>V + Vil {Vi — Vi) + ViflV ~ Vi) + b Vin{Vi ~ V n ) 

= (yiO + Vil + Vi2 + • • • + Vin)Vi - s/ii Vi - y i2 V 2 yinV l (6.23) 


or 

n n 

h ~ Vi ^2 yij — yij Vj 

3=0 3=1 


3 # i 


(6.24) 
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Vi Vi 

I Vil I 



FIGURE 6.7 

A lypical bus of the power system. 


The real and reactive power at bus i is 


P i +jQi = V i IT 


(6.25) 


Pi - jQi 


(6.26) 


Substituting for i* in (6.24) yields 


Pi- jQi _ „V' ^ T T 

Y* V Vij l^VijVj 
i 3 = O i=l 


(6.27) 


From the above relation, the mathematical formulation of the power flow 
problem results in a system of algebraic nonlinear equations which must be solved 
by iterative techniques. 


6.5 GAUSS-SEIDEL POWER FLOW SOLUTION 


In the power flow study, it is necessary to solve the set of nonlinear equations 
represented by (6.27) for two unknown variables at each node. In the Gauss-Seidel 
method (6.27) is solved for Vi, and the iterative sequence becomes 


Pf^-iOr* , r-. F (fc) 
■‘-.i*:) +L.y*j v j 


(6.28) 
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where yy shown in lowercase letters is the actual admittance in per unit. Pf ch 
and Qf ch are the net real and reactive powers expressed in per unit. In writing the 
KCL, current entering bus i was assumed positive. Thus, for buses where real and 
reactive powers are injected into the bus, such as generator buses, P? ch and Qf ch 
have positive values. For load buses where real and reactive powers are flowing 
away from the bus, P? ch and Q\ ch have negative values. If (6.27) is solved for P x 
and Qi, we have 

if +1) = 5 [vf jr J/y - jr yij v} k) ]} j ± i (6.29) 

3 = 0 3 - 1 

Q\ k+1) = -W <fc) [vf } jr, Vij - ± yy vf } ]} j /- * (6.30) 

j = 0 3 = 1 


The power flow equation is usually expressed in terms of the elements of 
the bus admittance matrix. Since the off-diagonal elements of the bus admittance 
matrix Yb us , shown by uppercase letters, are Yy = — yy, and the diagonal elements 
are Y ix = J2 Vij » (6.28) becomes 


pach _ 


iQi _ 


Vi 


(fc+i) _ 


V: 


Jk) 


E 




km 


(k) 


Yu 


(6.31) 


and 


P? +l) = K{V* W [V^Yu + jr Y tj v} k) ]} j ± i (6.32) 

i#* 

Q (k+1) = [v} k) Yii + jr YijYj k) }} 3 + i (6.33) 

j=i 


Yu includes the admittance to ground of line charging susceptance and any other 
fixed admittance to ground. In Section 6.7, a model is presented for transformers 
containing off-nominal ratio, which includes the effect of transformer tap setting. 

Since both components of voltage are specified for the slack bus, there are 
2 (n — 1) equations which must be solved by an iterative method. Under normal 
operating conditions, the voltage magnitude of buses are in the neighborhood of 
1 .0 per unit or close to the voltage magnitude of the slack bus. Voltage magnitude at 
load buses are somewhat lower than the slack bus value, depending on the reactive 
power demand, whereas the scheduled voltage at the generator buses are somewhat 
higher. Also, the jphase angle of the load buses are below the reference angle in 
accordance to the real power demand, whereas the phase angle of the generator 
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buses may be above the reference value depending on the amount of real power 
flowing into the bus. Thus, for the Gauss-Seidel method, an initial voltage estimate 
of 1.0 + jO.O for unknown voltages is satisfactory, and the converged solution 
correlates with the actual operating States. 

For P-Q buses, the real and reactive powers Pf ch and Qf ch are known. Start- 
ing with an initial estimate, (6.31) is solved for the real and imaginary components 
of voltage. For the voltage-controlled buses (P-V buses) where P? ch and |V*| are 
specified, lirst (6.33) is solved for Q\ k+1 \ and then is used in (6.31) to solve for 
V'/ fc+1) . However, since \Vi\ is specified, only the imaginary part of v} k+1 ^ is re- 
tained, and its real part is selected in order to satisfy 

( e ( fc + 1) )2 + (/ ( fc + D )2 = | y .|2 ( 6 . 34 ) 


or 

e\ k+1) = \J\Vi \ 2 — (// fe+1) ) 2 (6.35) 

where e| fc+1 ^ and f\ k+ ^ are the real and imaginary components of the voltage 
vf +1) in the iterative sequence. 

The rate of convergence is increased by applying an acceleration factor to the 
approximate solution obtained from each iteration. 

Vf +1) = V t ik) + a{V^ al - V} k) ) (6.36) 


where a is the acceleration factor. Its value depends upon the system. The range of 
1.3 to 1.7 is found to be satisfactory for typical systems. 

The updated voltages immediately replace the previous values in the solution 
of the subsequent equations. The process is continued until changes in the real and 
imaginary components of bus voltages between successive iterations are within a 
specified accuracy, i.e.. 


|ef +1 > 

l/f +1) 


-«n<« 

- 


(6.37) 


For the power mismatch to be reasonably small and acceptable, a very tight tol- 
erance must be specified on both components of the voltage. A voltage accuracy 
in the range of 0.00001 to 0.00005 pu is satisfactory. In practice, the method for 
determining the completion of a solution is based on an accuracy index set up on 
the power mismatch. The iteration continues until the magnitude of the largest ele- 
ment in the A P and A Q columns is less than the specified value. A typical power 
mismatch accuracy is 0.001 pu 

Once a solution is converged, the riet real and reactive powers at the slack bus 
are computed from (6.32) and (6.33). 
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6.6 LINE FLOWS AND LOSSES 

After the iterative solution of bus voltages, the next step is the computation of line 
flows and line losses. Consider the line connecting the two buses i and j in Figure 
6.8. The line current Iij, measured at bus i and defined positive in the direction 



FIGURE 6.8 

Transmission line model for calculating line flows. 
i —> j is given by 

hj — h + ho — Vij(Vi — Vj) + VioVi (6.38) 

Similarly, the line current Iji measured at bus j and defined positive in the direction 
j —y i is given by 

Iji — —h + Ijo = Vij{Vj — Vi) + VjoVj (6.39) 

The complex powers Sij from bus i to j and Sji from bus j to i are 

Sij = Vatj (6.40) 

Sji = Vjlji (6.41) 

The power loss in line z — j is the algebraic sum of the power flows determined 
from (6.40) and (6.41), i.e., 

S l ij — Sij + Sji (6.42) 

The power flow solution by the Gauss-Seidel method is demonstrated in the 
following two examples. 

Example 6.7 

Figure 6.9 shows the one-line diagram of a simple three-bus power system with 
generation at bus 1. The magnitude of voltage at bus 1 is adjusted to 1.05 per 
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unit. The scheduled loads at buses 2 and 3 are as marked on the diagram. Line 
impedances are marked in per unit on a 100-MVA base and the line charging sus- 
ceptances are neglected. 



138.6 45.2 

MW Mvar 


FIGURE 6.9 

One-line diagram of Example 6.7 (impedances in pu on 100-MVA base). 

(a) Using the Gauss-Seidel method, determine the phasor values of the voltage at 
the load buses 2 and 3 (P-Q buses) accurate to four decimal places. 

(b) Find the slack bus real and reactive power. 

(c) Determine the line flows and line losses. Construct a power flow diagram show- 
ing the direction of line flow. 


(a) Line impedances are converted to admittances 


y 12 = 


= 10 - j20 


" 0.02 + .70.04 

Similarly, j /13 = 10 — j'30 and 7/23 = 16 — j32. The admittances are marked on the 
network shown in Figure 6.10. 

At the P-Q buses, the complex loads expressed in per units are 

Sf 1 = _ ( 256 - 6 + J' 110 - 2 ) = -2.566 - jl.102 pu 
2 100 

5£* = -. ( 13 - . - 6 til 5 ; 2 ) = _ L386 _ jo.452 pu 
3 100 jr 

Since the actual admittances are readily available in Figure 6.10, for hånd calcu- 

lation, we use (6.28). Bus 1 is taken as reference bus (slack bus). Starting from 

an initial estimate of = 1.0 + j' 0.0 and = 1.0 + j' 0 . 0 . Vi and V 3 are 

computed from (6.28) as follows 

pj ch -m ch t /( q ) 


= -2.566- jl.102 pu 
= -1.386 -jO.452 pu 


= 


- + y 12 Vi + J/23V3 

V12 + 2/23 
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2/12 = 10 - j 20 


. 256.6 
MW 


2/13 = 10 - j 30 


2/23 = 16 — i 32 


110.2 

Mvar 


Slack Bus 3 -H Ly- 

Vi = 1.05/0° 

138.6 45.2 

MW Mvar 

FIGURE 6.10 

One-line diagram of Example 6.7 (admittances in pu on 100-MVA base). 

— 2.566+jl.l02 , /in „-onWi nc i i /i c 


_ + (10 ~ j 20) (1.05 + j 0) + (16 - j32)(1.0 + j 0) 


(26 - j 52) 


= 0.9825 - j'0.0310 


psch _Ar\sch /i\ 

J dm* + + tøVf > 

= LI 

2/13 + 2/23 

_ - 1 - 38 i 6 _+/° 452 + (io - j'30)(1.05 + jO) + (16 - j'32)(0.9825 - j0.0310) 

= (26 - j62) 

= 1.0011 -jO.0353 

For the second iteration we have 

(2 ) . i ggo + ( 10 - J 20) (1.05 + jO) + (16 - j32)(1.0011 - j0.0353) 
2 (26 - jr'52) 

= 0.9816 - j'0.0520 


(2) looi 8 Æ^53 + (10 - j"30)(1.05 + jO) + (16 - i32)(0.9816 - ;0.052) 
3 (26 - j62) 

= 1.0008- jO. 0459 

The process is continued and a solution is converged with an accuracy of 5 x 10 -5 
per unit in seven iterations as given below. 

V 2 (3) = 0.9808 - jO.0578 F 3 (3) = 1.0004 - j0.0488 
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F 2 (4) = 0.9803 - j0.0594 
K 2 (5) = 0.9801 - jO.0598 
V 2 (6) = 0.9801 - jO.0599 
vj 7) = 0.9800 - 7O.O6OO 

The final solution is 


V 3 (4) = 1.0002 - jO.0497 
V 3 (5) = 1.0001 - 7 0.0499 
V 3 (6) = 1.0000 - jO.0500 
V 3 (7) = 1.0000 - jO.0500 


V 2 = 0.9800 - 7O.O6OO = 0.98183Z -3.5035° pu 
V 3 = 1.0000 - jO.0500 = 1.00125Z -2.8624° pu 


(b) With the knowledge of all bus voltages, the slack bus power is obtained from 
(6.27) 

Pi - jQi = V* [Vi (yi2 + 2/13) - (2/12V2 + 2/13V3)] 

= 1.05(1.05(20 - j50) - (10 - 7 20) (0.98 - j.06) - 
(10 — j30)(1.0 — j'0.05)] 

= 4.095 -jl.890 


or the slack bus real and reactive powers are Pi = 4.095 pu = 409.5 MW and 
Qi = 1.890 pu = 189 Mvar. 


(c) To find the line flows, first the line currents are computed. With line charging 
capacitors neglected, the line currents are 


ha = Vi 2 (Vi - V 2 ) = (10 - 7 20) [(1.05 + 7O) - (0.98 - j0.06)] = 1.9 - j0.8 
I 21 = —I \2 — —1-9 + 7O.8 


hs = Vis(Vi - V 3 ) = (10 - 730) [(1.05 + j0) - (1.0 - 70.05)] = 2.0 - 7I.O 

1 31 = -ii3 = —2.0 + 7I.O 

las = m(Va - V 3 ) = (16 - j32)[(0.98 - 7O.O6) - (1 - j0.05)] = -.64 + jA 8 

1 32 = — Z23 = 0.64 — 7'0.48 

The line flows are 


S12 = VJu = (1.05 + 7'0.0)(1.9 + 7O.8) = 1.995 + j0.84 pu 
= 199.5 MW + 784.0 Mvar 

S 2 i = V 2 I 21 = (0.98 - 7'0.06)(— 1.9 - 7O.8) = -1.91 - j0.67 pu 
= -191.0 MW- 767.O Mvar 
Sis = Vil { 3 = (1.05 + j0.0)(2.0 + 7I O) = 2.1 + 71.05 pu 
= 210.0 MW + 7105.0 Mvar 
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5 3 1 = V 3 In = (1.0 - j0.05)(-2.0 - j 1.0) = -2.05 - i0.90 pu 

= —205.0 MW — j90.0 Mvar 

S 2 z = V 2/23 = (0-98 - j0.06)(— 0.656 + j0.48) = -0.656 - jO.432 pu 
= -65.6 MW - j43.2 Mvar 

5 32 = V 3 I 32 = (1-0 - j0.05)(0.64 + j'0.48) = 0.664 + j0.448 pu 

= 66.4 MW + j44.8 Mvar 

and the line losses are 

Sl 12 = 5x2 + 5 2 i = 8.5 MW + jl7.0 Mvar 

Sl 13 = 5 x 3 + 5 3 x = 5.0 MW + J15.0 Mvar 

Sl 23 = ^23 + S 32 — 0.8 MW + j 1 .60 Mvar 

The power flow diagram is shown in Figure 6.11, where real power direction is 

indicated by — > and the reactive power direction is indicated by 1 — >. The values 
within parentheses are the real and reactive losses in the line. 



FIGURE 6.11 

Power flow diagram of Example 6.7 (powers in MW and Mvar). 


Example 6.8 

Figure 6.12 shows the one-line diagram of a simple three-bus power system with 
generators at buses 1 and 3. The magnitude of voltage at bus 1 is adjusted to 1.05 
pu. Voltage magnitude at bus 3 is fixed at 1.04 pu with a real power generation 
of 200 MW. A load consisting of 400 MW and 250 Mvar is taken from bus 2. 
Line impedances are marked in per unit on a 100 MVA base, and the line charging 
susceptances are neglected. Obtain the power flow solution by the Gauss-Seidel 
method including line flows and line losses. 
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FIGURE 6.12 

One-line diagram of Example 6.8 (impedances in pu on 100-MVA base). 
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The value of Q' 3 ’ is used as Qf h for the computation of voltage at bus 3. The 
complex voltage at bus 3, denoted by V^\ is calculated 

r>9ch__Ar\sch /n\ 

Jft - + yuV, + m,V. 

yUJ v 3 

c3 2/13 + 2/23 

_ 2 i°q 4 ^q - + (10 - J30)(1.05 + jO) + (16 - j'32) (0.97462 - j0.042307) 


(26 - j62) 


= 1.03783 - jO. 005170 


Since IV 3 I is held constant at 1.04 pu, only the imaginary part of is retained, 
i.e, = —0.005170, and its real part is obtained from 

eg 1} = 7(1.04)2 - (0.005170) 2 = 1.039987 


V 3 (1) = 1.039987 — y'0.005170 
For the second iteration, we have 
+ yi2 v t + 

Vj ’ = 12 

2/12 + 2/23 

_ : 97 ^62+^0423b7 + ( 10 - j 20)(1.05) + (16 - j32)(1.039987 + j0.005170) 


(26 - j52) 


0.971057 - jO.043432 


Q? ] = -nvi (1) [y 3 ( 1 ) ( 2 /l 3 + 2/23) - 2 / 13 Vi - 2 / 23 V 2 (2) ]} 

= -9{(1.039987 + j0.005170)[(1.039987 - j0.005170)(26 - j62) - 
(10 - j 30) (1.05 + jO) - (16 - j32)(0. 971057 - jO.043432)]} 

= 1.38796 


psch __Af)Sch (<y\ 

2/13 Vi + 2/23 v 2 (2) 

_ V 3 

2/13 + 2/23 

llSSlT + ( 10 - 7’30)(1.05) + (16 - j32)(. 971057 - j.043432) 

(26 — j62) 


= 1.03908 - j'0.00730 
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Since | V 3 1 is held constant at 1.04 pu, only the imaginary part of is retained, 

( 2 ) 

i.e, fl = —0.00730, and its real part is obtained from 

e!j 2) = \j (1.04) 2 - (0.00730) 2 = 1.039974 
or 

y 3 (2) = 1.039974 - jO. 00730 


The process is continued and a solution is converged with an accuracy of 5 x 10 5 
pu in seven iterations as given below. 


CO 

= 0.97073 - y'0.04479 

Q ( 3 3) = 1.42904 

II 

'S? 

1.03996 - jf0.00833 


= 0.97065 - jO.04533 

Q ( 3 4) = 1.44833 

II 

"S” 

1.03996 -j 0.00873 

vf* 

= 0.97062 - jO.04555 

QP = 1.45621 

V 3 ,s> = 

1.03996 -j 0.00893 

Vf> 

= 0.97061 - jO.04565 

Q^ 6) = 1.45947 

vf - 

1.03996 - j0.00900 

v#> 

= 0.97061 - jO.04569 

QP = 1.46082 

n' 71 = 

1.03996 - j0.00903 

The final solution is 





V 2 = 

0.97168Z— 2.6948° 

pu 



S 3 = 

2.0 + ^1.4617 pu 




^3 = 

1.04Z — .498° pu 




Sl = 

2.1842 + j 1.4085 

pu 



Line flows and line losses are computed as in Example 6.7, and the results ex- 
pressed in MW and Mvar are 

5i2 = 179.36 + jll8.734 S 2 i = -170.97 - jl01.947 S L12 = 8.39 + y'16.79 

S = 39.06 + j 22. 118 5 3 i = -38.88 - j 21.569 S L 13 = 0.18 + j0.548 

S 2 3 = -229.03 - jl48.05 S 32 = 238.88 + jl67.746 S L2 3 = 9.85 + jl9M 

The power flow diagram is shown in Figure 6.13, where real power direction 
is indicated by — > and the reactive power direction is indicated by 1 — >. The values 
within parentheses are the real and reactive losses in the line. 
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FIGURE 6.13 

Power flow diagram of Example 6.8 (powers in MW and Mvar). 


6.7 TAP CHANGING TRANSFORMERS 

In Section 2.6 it was shown that the flow of real power along a transmission line is 
determined by the angle difference of the terminal voltages, and the flow of reactive 
power is determined mainly by the magnitude difference of terminal voltages. Real 
and reactive powers can be controlled by use of tap changing transformers and 
regulating transformers. 

In a tap changing transformer, when the ratio is at the nominal value, the 
transformer is represented by a series admittance yt in per unit. With off-nominal 
ratio, the per unit admittance is different from both sides of the transformer, and the 
admittance must be modified to include the effect of the off-nominal ratio. Consider 
a transformer with admittance y t in series with an ideal transformer representing 
the off-nominal tap ratio l:o as shown in Figure 6.14. yt is the admittance in per 
unit based on the nominal tum ratio and a is the per unit off-nominal tap position 
allowing for small adjustment in voltage of usually ±10 percent. In the case of 
phase shifting transformers, a is a complex number. Consider a fictitious bus x 
between the tum ratio and admittance of the transformer. Since the complex power 
on either side of the ideal transformer is the same, it follows that if the voltage goes 
through a positive phase angle shift, the current will go through a negative phase 
angle shift. Thus, for the assumed direction of currents, we have 

^ = ~Vj (6.43) 

a J 

h = —a*Ij (6.44) 

The current I t is given by 
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Vj 


F1GURE 6.14 

Transformer with tap setting ratio a:l 


Substituting for V x , we have 


li = ytVi - y -Vj 
a 


(6.45) 


Also, from (6.44) we have 



substituting for li from (6.45) we have 



+ -iLy. 

+ lal 2 J 


writing (6.45) and (6.46) in matrix form results in 





(6.46) 


(6.47) 


For the case when a is real, the 7r model shown in Figure 6.15 represents the ad- 
mittance matrix in (6.47). In the 7r model, the left side corresponds to the non-tap 
side and the right side corresponds to the tap side of the transformer. 



FIGURE 6.15 

Equivalent Circuit for a tap changing transformer. 
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6.8 POWER FLOW PROGRAMS 

Several computer programs have been developed for the power flow solution of 
practical systems. Each method of solution consists of four programs. The pro- 
gram for the Gauss-Seidel method is Ifgauss, which is preceded by lfybus, and is 
followed by busout and lineflow. Programs lfybus, busout, and lineflow are de- 
signed to be used with two more power flow programs. These are Ifnewton for 
the Newton-Raphson method and decouple for the fast decoupled method. The 
following is a brief description of the programs used in the Gauss-Seidel method. 


lfybus This program requires the line and transformer parameters and transformer 
tap settings specified in the input file named linedata. It converts impedances 
to admittances and obtains the bus admittance matrix. The program is de- 
signed to handle parallel lines. 


Ifgauss This program obtains the power flow solution by the Gauss-Seidel method 
and requires the files named busdata and linedata. It is designed for the di- 
rect use of load and generation in MW and Mvar, bus voltages in per unit, 
and angle in degrees. Loads and generation are converted to per unit quanti- 
ties on the base MVA selected. A provision is made to maintain the generator 
reactive power of the voltage-controlled buses within their specified limits. 
The violation of reactive power limit may occur if the specified voltage is 
either too high or too low. After a few iterations (10 f, ‘ iteration in the Gauss 
method), the var calculated at the generator buses are examined. If a limit is 
reached, the voltage magnitude is adjusted in steps of 0.5 percent up to ±5 
percent to bring the var demand within the specified limits. 


busout This program produces the bus output result in a tabulated form. The bus 
output result includes the voltage magnitude and angle, real and reactive 
power of generators and loads, and the shunt capacitor/reactor Mvar. Total 
generation and total load are also included as outlined in the sample case. 


lineflow This program prepares the line output data. It is designed to display the 
active and reactive power flow entering the line terminals and line losses as 
well as the net power at each bus. Also included are the total real and reactive 
losses in the system. The output of this portion is also shown in the sample 
case. 
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6.9 DATA PREPARATION 

In order to perform a power flow analysis by the Gauss-Seidel method in the MAT- 
LAB environment, the following variables must be defined: power system base 
MVA, power mismatch accuracy, acceleration factor, and maximum number of it- 
erations. The name (in lowercase letters) reserved for these variables are basem va, 
accuracy, accel, and maxiter, respectively. Typical values are as follows: 

basemva = 100; accuracy = 0.001; 
accel = 1.6; maxiter = 80; 

The initial step in the preparation of input file is the numbering of each bus. Buses 
are numbered sequentially. Although the numbers are sequentially assigned, the 
buses need not be entered in sequence. In addition, the following data files are re- 
quired. 

BUS DATA FILE - busdata The format for the bus entry is chosen to facili- 
tate the required data for each bus in a single row. The information required must be 
included in a matrix called busdata. Column 1 is the bus number. Column 2 con- 
tains the bus code. Columns 3 and 4 are voltage magnitude in per unit and phase 
angle in degrees. Columns 5 and 6 are load MW and Mvar. Column 7 through 10 
are MW, Mvar, minimum Mvar and maximum Mvar of generation, in that order. 
The last column is the injected Mvar of shunt capacitors. The bus code entered in 
column 2 is used for identifying load, voltage-controlled, and slack buses as out- 
lined below: 

1 This code is used for the slack bus. The only necessary information for this bus 

is the voltage magnitude and its phase angle. 

0 This code is used for load buses. The loads are entered positive in megawatts 
and megavars. For this bus, initial voltage estimate must be specified. This is 
usually 1 and 0 for voltage magnitude and phase angle, respectively. If volt- 
age magnitude and phase angle for this type of bus are specified, they will 
be taken as the initial starting voltage for that bus instead of a flat start of 1 
and 0. 

2 This code is used for the voltage-controlled buses. For this bus, voltage magni- 

tude, real power generation in megawatts, and the minimum and maximum 
limits of the megavar demand must be specified. 

LINE DATA FILE - linedata Lines are identified by the node-pair method. The 
information required must be included in a matrix called linedata. Columns 1 and 
2 are the line bus numbers. Columns 3 through 5 contain the line resistance, reac- 
tance, and one-half of the total line charging susceptance in per unit on the specified 
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MVA base. The last column is for the transformer tap setting; for lines, 1 must be 
entered in this column. The lines may be entered in any sequence or order with 
the only restriction being that if the entry is a transformer, the left bus number is 
assumed to be the tap side of the transformer. 

The IEEE 30 bus system is used to demonstrate the data preparation and the 
use of the power flow programs by the Gauss-Seidel method. 

Example 6.9 

Figure 6.16 is part of the American Electric Power Service Corporation network 
which is being made available to the electric Utility industry as a standard test case 
for evaluating various analytical methods and computer programs for the solution 
of power system problems. Use the lfgauss program to obtain the power solution 
by the Gauss-Seidel method. Bus 1 is taken as the slack bus with its voltage ad- 
justed to 1.06Z0 0 pu. The data for the voltage-controlled buses is 


Regulated Bus Data 

Bus 

Voltage 

Min. Mvar 

Max. Mvar 

No. 

Magnitude 

Capacity 

Capacity 

2 

1.043 

-40 

HHE ■ 

5 

1.010 

-40 


8 

1.010 

-10 


11 

1.082 

-6 


13 

1.071 

-6 



Transformer tap setting are given in the table below. The left bus number is as- 
sumed to be the tap side of the transformer. 



The data for the injected Q due to shunt capacitors is 


Injected Q due to Capacitors 

Bus No. Mvar 

_ _ - 

24 4.3 
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1 



FIGURE 6.16 

30-bus IEEE sample system. 
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Generation and loads are as given in the data prepared for use in the MATLAB 
environment in the matrix defined as busdata. Code 0, code 1, and code 2 are used 
for the load buses, the slack bus and the voltage-controlled buses, respectively. 
Values for basemva, accuracy, accel and maxiter must be specified. Line data are 
as given in the matrix called linedata. The last column of this data must contain 1 
for lines, or the tap setting values for transformers with off-nominal tum ratio. The 
control commands required are Ifybus, Ifgauss and lineflow. A diary command 
may be used to save the output to the specified file name. The power flow data and 
the commands required are as follows. 

clear ’/, clears all variables from workspace . 

basemva = 100; accuracy = 0.001; accel = 1.8; maxiter = 100; 

*/. IEEE 30-BUS TEST SYSTEM (American Electric Power) 

7, Bus Bus Voltage Angle — Load — Generator Injected 

7, No code Mag. Degree MW Mvar MW Mvar Qmin Qmax Mvar 


bus dat a= [1 

1 

1.06 

0 

0.0 

0.0 

0.0 

0.0 

0 

0 

0 

2 

2 

1.043 

0 

21.70 

12.7 

40.0 

0.0 

-40 

50 

0 

3 

0 

1.0 

0 

2.4 

1.2 

0.0 

0.0 

0 

0 

0 

4 

0 

1.06 

0 

7.6 

1.6 

0.0 

0.0 

0 

0 

0 

5 

2 

1.01 

0 

94.2 

19.0 

0.0 

0.0 

-40 

40 

0 

6 

0 

1.0 

0 

0.0 

0.0 

0.0 

0.0 

0 

0 

0 

7 

0 

1.0 

0 

22.8 

10.9 

0.0 

0.0 

0 

0 

0 

8 

2 

1.01 

0 

30.0 

30.0 

0.0 

0.0 

-10 

40 

0 

9 

0 

1.0 

0 

0.0 

0.0 

0.0 

0.0 

0 

0 

0 

10 

0 

1.0 

0 

5.8 

2.0 

0.0 

0.0 

0 

0 

19 

11 

2 

1.082 

0 

0.0 

0.0 

0.0 

0.0 

-6 

24 

0 

12 

0 

1.0 

0 

11.2 

7.5 

0 

0 

0 

0 

0 

13 

2 

1.071 

0 

0.0 

0.0 

0 

0 

-6 

24 

0 

14 

0 

1.0 

0 

6.2 

1.6 

0 

0 

0 

0 

0 

15 

0 

1.0 

0 

8.2 

2.5 

0 

0 

0 

0 

0 

16 

0 

1.0 

0 

3.5 

1.8 

0 

0 

0 

0 

0 

17 

0 

1.0 

0 

9.0 

5.8 

0 

0 

0 

0 

0 

18 

0 

1.0 

0 

3.2 

0.9 

0 

0 

0 

0 

0 

19 

0 

1.0 

0 

9.5 

3.4 

0 

0 

0 

0 

0 

20 

0 

1.0 

0 

2.2 

0.7 

0 

0 

0 

0 

0 

21 

0 

1.0 

0 

17.5 

11.2 

0 

0 

0 

0 

0 

22 

0 

1.0 

0 

0.0 

0.0 

0 

0 

0 

0 

0 

23 

0 

1.0 

0 

3.2 

1.6 

0 

0 

0 

0 

0 

24 

0 

1.0 

0 

8.7 

6.7 

0 

0 

0 

0 

4.3 

25 

0 

1.0 

0 

0.0 

0.0 

0 

0 

0 

0 

0 

26 

0 

1.0 

0 

3.5 

2.3 

0 

0 

0 

0 

0 

27 

0 

1.0 

0 

0.0 

0.0 

0 

0 

0 

0 

0 

28 

0 

1.0 

0 

0.0 

0.0 

0 

0 

0 

0 

0 

29 

0 

1.0 

0 

2.4 

0.9 

0 

0 

0 

0 

0 

30 

0 

1.0 

0 

10.6 

1.9 

0 

0 

0 

0 

0]; 
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'/, Line Data 
% 


y, Bus 

bus 

R 

X 

1/2 B 1 for Line code or 

7, ni 

nr 

pu 

pu 

pu tap setting value 

linedata=[l 

2 

0.0192 

0 . 0575 

0.02640 

1 

1 

3 

0 . 0452 

0.1852 

0.02040 

1 

2 

4 

0.0570 

0.1737 

0.01840 

1 

3 

4 

0.0132 

0.0379 

0 . 00420 

1 

2 

5 

0.0472 

0.1983 

0.02090 

1 

2 

6 

0.0581 

0.1763 

0.01870 

1 

4 

6 

0.0119 

0.0414 

0.00450 

1 

5 

7 

0.0460 

0.1160 

0.01020 

1 

6 

7 

0.0267 

0.0820 

0.00850 

1 

6 

8 

0.0120 

0 . 0420 

0.00450 

1 

6 

9 

0.0 

0.2080 

0.0 

0.978 

6 

10 

0.0 

0.5560 

0.0 

0.969 

9 

11 

0.0 

0.2080 

0.0 

1 

9 

10 

0.0 

0.1100 

0.0 

1 

4 

12 

0.0 

0.2560 

0.0 

0.932 

12 

13 

0.0 

0.1400 

0.0 

1 

12 

14 

0.1231 

0.2559 

0.0 

1 

12 

15 

0.0662 

0.1304 

0.0 

1 

12 

16 

0.0945 

0.1987 

0.0 

1 

14 

15 

0.2210 

0.1997 

0.0 

1 

16 

17 

0.0824 

0.1923 

0.0 

1 

15 

18 

0.1073 

0.2185 

0.0 

1 

18 

19 

0.0639 

0.1292 

0.0 

1 

19 

20 

0 . 0340 

0.0680 

0.0 

1 

10 

20 

0.0936 

0.2090 

0.0 

1 

10 

17 

0 . 0324 

0.0845 

0.0 

1 

10 

21 

0.0348 

0 . 0749 

0.0 

1 

10 

22 

0.0727 

0.1499 

0.0 

1 

21 

22 

0.0116 

0.0236 

0.0 

1 

15 

23 

0.1000 

0.2020 

0.0 

1 

22 

24 

0.1150 

0.1790 

0.0 

1 

23 

24 

0.1320 

0.2700 

0.0 

1 

24 

25 

0.1885 

0.3292 

0.0 

1 

25 

26 

0 . 2544 

0.3800 

0.0 

1 

25 

27 

0.1093 

0.2087 

0.0 

1 

28 

27 

0.0000 

0.3960 

0.0 

0.968 

27 

29 

0.2198 

0.4153 

0.0 

1 

27 

30 

0.3202 

0.6027 

0.0 

1 

29 

30 

0.2399 

0.4533 

0.0 

1 

8 

28 

0.0636 

0 . 2000 

0.0214 

1 

6 

28 

0.0169 

0.0599 

0.065 

i]; 
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lfybus 
lf gauss 
busout 
linef low 


% Forms the bus admittance matrix 
7. Power flow solution by Gauss-Seidel method 
7, Prints the power flow solution on the screen 
7. Computes and displays the line flow and losses 


The Ifgauss, busout and the lineflow produce the following tabulated results. 


Power Flow Solution by Gauss-Seidel Method 
Maximum Power mismatch = 0.000951884 
No. of iterations = 34 


Bus 

Voltage 

s Angle 

Load 

— Generation — Injected 

No. 

Mag. 

Degree 

MW 

Mvar 

MW 

Mvar 

Mvar 

1 

1.060 

0.000 

0.000 

0.000 

260.950 

-17.010 

0.00 

2 

1.043 

-5.496 

21.700 

12.700 

40 . 000 

48.826 

0.00 

3 

1.022 

-8 . 002 

2.400 

1.200 

0.000 

0.000 

0.00 

4 

1.013 

-9.659 

7.600 

1.600 

0.000 

0.000 

0.00 

5 

1.010 

-14.380 

94.200 

19.000 

0.000 

35.995 

0.00 

6 

1.012 

-11.396 

0.000 

0.000 

0.000 

0.000 

0.00 

7 

1.003 

-13.149 

22.800 

10.900 

0.000 

0.000 

0.00 

8 

1.010 

-12.114 

30.000 

30.000 

0.000 

30.759 

0.00 

9 

1.051 

-14.432 

0.000 

0.000 

0.000 

0.000 

0.00 

10 

1.044 

-16.024 

5.800 

2.000 

0.000 

0.000 

19.00 

11 

1.082 

-14.432 

0.000 

0.000 

0.000 

16.113 

0.00 

12 

1.057 

-15.301 

11.200 

7.500 

0.000 

0.000 

0.00 

13 

1.071 

-15.300 

0.000 

0.000 

0.000 

10.406 

0.00 

14 

1.043 

-16.190 

6.200 

1.600 

0.000 

0.000 

0.00 

15 

1.038 

-16.276 

8.200 

2.500 

0.000 

0.000 

0.00 

16 

1.045 

-15.879 

3.500 

1.800 

0.000 

0.000 

0.00 

17 

1.039 

-16.187 

9.000 

5.800 

0.000 

0.000 

0.00 

18 

1.028 

-16.881 

3.200 

0.900 

0.000 

0.000 

0.00 

19 

1.025 

-17.049 

9.500 

3.400 

0.000 

0.000 

0.00 

20 

1.029 

-16.851 

2.200 

0.700 

0.000 

0.000 

0.00 

21 

1.032 

-16.468 

17.500 

11.200 

0.000 

0.000 

0.00 

22 

1.033 

-16.455 

0.000 

0.000 

0.000 

0.000 

0.00 

23 

1.027 

-16.660 

3.200 

1.600 

0.000 

0.000 

0.00 

24 

1.022 

-16.829 

8.700 

6.700 

0.000 

0.000 

4.30 

25 

1.019 

-16.423 

0.000 

0.000 

0.000 

0.000 

0.00 

26 

1.001 

-16.835 

3.500 

2.300 

0.000 

0.000 

0.00 

27 

1.026 

-15.913 

0.000 

0.000 

0.000 

0.000 

0.00 

28 

1.011 

-12.056 

0.000 

0.000 

0.000 

0.000 

0.00 
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29 1 . 006 

-17.133 

2.400 

0.900 

0.000 

0.000 

0.00 

30 0.994 

-18.016 

10.600 

1.900 

0.000 

0.000 

0.00 

Total 


283.400 

126.200 

300.950 

125.089 

23.30 


Line Flow and Losses 


— Line — 

Power at 

bus & line flow 

— Line loss — Transformer 

from 

to 

MW 

Mvar 

MVA 

MW 

Mvar 

tap 

1 


260.950 

-17.010 

261.504 





2 

177.743 

-22.140 

179.117 

5.461 

10.517 



3 

83.197 

5.125 

83.354 

2.807 

7.079 


2 


18.300 

36.126 

40.497 





1 

-172.282 

32.657 

175.350 

5.461 

10.517 



4 

45.702 

2.720 

45.783 

1.106 

-0.519 



5 

82.990 

1.704 

83.008 

2.995 

8.178 



6 

61.905 

-0.966 

61.913 

2.047 

2.263 


3 


-2.400 

-1.200 

2.683 





1 

-80.390 

1.954 

80.414 

2.807 

7.079 



4 

78 . 034 

-3.087 

78.095 

0.771 

1.345 


4 


-7.600 

-1.600 

7.767 





2 

-44.596 

-3.239 

44.713 

1.106 

-0.519 



3 

-77.263 

4.432 

77.390 

0.771 

1.345 



6 

70.132 

-17.624 

72.313 

0.605 

1.181 



12 

44.131 

14.627 

46.492 

0.000 

4.686 

0.932 

5 


-94.200 

16.995 

95.721 





2 

-79.995 

6.474 

80.256 

2.995 

8.178 



7 

-14.210 

10.467 

17.649 

0.151 

-1.687 


6 


0.000 

0.000 

0.000 





2 

-59.858 

3.229 

59.945 

2.047 

2.263 



4 

-69.527 

18.805 

72.026 

0.605 

1.181 



7 

37.537 

-1.915 

37.586 

0.368 

-0.598 



8 

29.534 

-3.712 

29.766 

0.103 

-0.558 



9 

27.687 

-7.318 

28.638 

0.000 

1.593 

0.978 


10 

15.828 

0.656 

15.842 

0.000 

1.279 

0.969 


28 

18.840 

-9.575 

21.134 

0.060 

-13.085 


7 


-22.800 

-10.900 

25.272 





5 

14.361 

-12.154 

18.814 

0.151 

-1.687 
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6 

- 37 . 

.170 

1 . 

,317 

37 . 

.193 

0 . 

,368 

- 0 . 

,598 

8 


- 30 . 

.000 

0 . 

.759 

30 . 

.010 






6 

- 29 , 

.431 

3 , 

.154 

29 . 

.599 

0 . 

,103 

- 0 . 

.558 


28 

- 0 . 

.570 

- 2 . 

.366 

2 

.433 

0 . 

,000 

- 4 . 

368 

9 


0 , 

.000 

0 . 

.000 

0 , 

.000 






6 

- 27 . 

.687 

8 . 

.911 

29 , 

.086 

0 , 

.000 

1 , 

,593 


11 

0 . 

.003 

- 15 , 

.653 

15 . 

.653 

- 0 . 

.000 

0 . 

.461 


10 

27 . 

.731 

6 . 

.747 

28 . 

.540 

0 . 

.000 

0 . 

.811 

10 


- 5 . 

.800 

17 , 

.000 

17 

.962 






6 

- 15 . 

.828 

0 . 

.623 

15 . 

.840 

0 . 

,000 

1 . 

.279 


9 

- 27 , 

.731 

- 5 , 

.936 

28 . 

.359 

0 . 

,000 

0 . 

,811 


20 

9 . 

.018 

3 , 

.569 

9 . 

.698 

0 , 

,081 

0 . 

,180 


17 

5 . 

.347 

4 . 

.393 

6 . 

.920 

0 . 

.014 

0 , 

.037 


21 

15 . 

.723 

9 . 

.846 

18 , 

.551 

0 . 

,110 

0 . 

.236 


22 

7 . 

.582 

4 . 

.487 

8 

.811 

0 , 

.052 

0 . 

, 107 

11 


0 . 

.000 

16 , 

.113 

16 , 

.113 






9 

- 0 , 

.003 

16 , 

.114 

16 , 

.114 

- 0 , 

.000 

0 , 

.461 

12 


- 11 . 

.200 

- 7 . 

.500 

13 . 

.479 






4 

- 44 , 

.131 

- 9 , 

.941 

45 . 

.237 

0 , 

.000 

4 . 

.686 


13 

- 0 , 

.021 

- 10 , 

.274 

10 , 

.274 

0 . 

,000 

0 . 

.132 


14 

7 . 

.852 

2 . 

.428 

8 . 

.219 

0 . 

.074 

0 , 

.155 


15 

17 . 

.852 

6 . 

.968 

19 . 

.164 

0 . 

.217 

0 , 

.428 


16 

7 . 

.206 

3 , 

.370 

7 , 

.955 

0 . 

.053 

0 . 

.112 

13 


0 . 

,000 

10 , 

.406 

10 , 

.406 






12 

0 . 

.021 

10 , 

.406 

10 , 

.406 

0 , 

.000 

0 . 

.132 

14 


- 6 . 

.200 

- 1 . 

.600 

6 , 

.403 






12 

- 7 . 

.778 

- 2 . 

.273 

8 . 

.103 

0 . 

.074 

0 , 

.155 


15 

1 . 

,592 

0 , 

.708 

1 , 

.742 

0 . 

.006 

0 . 

.006 

15 


- 8 . 

.200 

- 2 , 

.500 

8 . 

.573 






12 

- 17 . 

,634 

- 6 , 

.540 

18 , 

.808 

0 . 

.217 

0 . 

.428 


14 

- 1 , 

.586 

- 0 . 

.702 

1 , 

.734 

0 . 

.006 

0 

.006 


18 

6 . 

,009 

1 , 

.741 

6 , 

.256 

0 . 

.039 

0 , 

.079 


23 

5 . 

.004 

2 . 

.963 

5 . 

.815 

0 . 

.031 

0 

.063 

16 


- 3 . 

.500 

- 1 , 

.800 

3 . 

.936 






12 

- 7 , 

.152 

- 3 . 

.257 

7 . 

.859 

0 . 

.053 

0 

.112 


17 

3 . 

,658 

1 . 

.440 

3 . 

.931 

0 . 

.012 

0 

.027 
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17 


- 9.000 

- 5.800 

10.707 




16 

- 3.646 

- 1.413 

3.910 

0.012 

0.027 


10 

- 5.332 

- 4.355 

6.885 

0.014 

0.037 

18 


- 3.200 

- 0.900 

3.324 




15 

- 5.970 

- 1.661 

6.197 

0.039 

0.079 


19 

2.779 

0.787 

2.888 

0.005 

0.010 

19 


- 9.500 

- 3.400 

10.090 




18 

- 2.774 

- 0.777 

2.881 

0.005 

0.010 


20 

- 6.703 

- 2.675 

7.217 

0.017 

0.034 

20 


- 2.200 

- 0.700 

2.309 




19 

6.720 

2.709 

7.245 

0.017 

0.034 


10 

- 8.937 

- 3.389 

9.558 

0.081 

0.180 

21 


- 17.500 

- 11.200 

20 . 777 




10 

- 15.613 

- 9.609 

18.333 

0.110 

0.236 


22 

- 1.849 

- 1.627 

2.463 

0.001 

0.001 

22 


0.000 

0.000 

0.000 




10 

- 7.531 

- 4.380 

8.712 

0.052 

0.107 


21 

1.850 

1.628 

2.464 

0.001 

0.001 


24 

5.643 

2.795 

6.297 

0.043 

0.067 

23 


- 3.200 

- 1.600 

3.578 




15 

- 4.972 

- 2.900 

5.756 

0.031 

0.063 


24 

1.771 

1.282 

2.186 

0.006 

0.012 

24 


- 8.700 

- 2.400 

9.025 




22 

- 5.601 

- 2.728 

6.230 

0.043 

0.067 


23 

- 1.765 

- 1.270 

2.174 

0.006 

0.012 


25 

- 1.322 

1.604 

2.079 

0.008 

0.014 

25 


0.000 

0.000 

0.000 




24 

1.330 

- 1.590 

2.073 

0.008 

0.014 


26 

3.520 

2.372 

4.244 

0.044 

0.066 


27 

- 4.866 

- 0.786 

4.929 

0.026 

0.049 

26 


- 3.500 

- 2.300 

4.188 




25 

- 3.476 

- 2.306 

4.171 

0.044 

0.066 

27 


0.000 

0.000 

0.000 




25 

4.892 

0.835 

4.963 

0.026 

0.049 
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28 

-18.192 

-4.152 


29 

6.178 

1.675 


30 

7.093 

1.663 

28 


0.000 

0.000 


27 

18.192 

5.463 


8 

0.570 

-2.003 


6 

-18.780 

-3.510 

29 


-2.400 

-0.900 


27 

-6.093 

-1.513 


30 

3.716 

0.601 

30 


-10.600 

-1.900 


27 

-6.932 

-1.359 


29 

-3.683 

-0.537 


Total loss 


18.660 

- 0.000 

1.310 

6.401 

0.086 

0.162 

7.286 

0.162 

0.304 

0.000 

18.994 

- 0.000 

1.310 

2.082 

0.000 

-4.368 

19.106 

0.060 

-13.085 

2.563 

6.278 

0.086 

0.162 

3.764 

0.034 

0.063 

10.769 

7.064 

0.162 

0.304 

3.722 

0.034 

0.063 


17.594 

22.233 


0.968 


6.10 NEWTON-RAPHSON POWER FLOW SOLUTION 

Because of its quadratic convergence, Newton’s method is mathematically superior 
to the Gauss-Seidel method and is less prone to divergence with ill-conditioned 
problems. For large power systems, the Newton-Raphson method is found to be 
more efficient and practical. The number of iterations required to obtain a solution 
is independent of the system size, but more functional evaluations are required at 
each iteration. Since in the power flow problem real power and voltage magnitude 
are specified for the voltage-controlled buses, the power flow equation is formu- 
lated in polar form. For the typical bus of the power system shown in Figure 6.7, 
the current entering bus i is given by (6.24). This equation can be rewritten in terms 
of the bus admittance matrix as 


= ( 6 - 4 8 ) 
3 = 1 

In the above equation, j includes bus i. Expressing this equation in polar form, we 
have 

h = f2\ Y vm^iJ + S j ( 6 - 49 ) 

j=i 


The complex power at bus i is 


Pi - jQi = vth 


(6.50) 
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Substituting from (6.49) for li in (6.50), 

Pi ~ jQi = \Vi\Z-6i £ \YyWVj\lOij + Sj (6.51) 

3=1 

Separating the real and imaginary parts, 

Pi = £ cos (% - 8i + Sj) (6.52) 

3=1 

Qi = - £ MlWilsin (% - ft + ft) (6.53) 

3=1 

Equations (6.52) and (6.53) constitute a set of nonlinear algebraic equations in 
terms of the independent variables, voltage magnitude in per unit, and phase angle 
in radians. We have two equations for each load bus, given by (6.52) and (6.53), and 
one equation for each voltage-controlled bus, given by (6.52). Expanding (6.52) 
and (6.53) in Taylor’s series about the initial estimate and neglecting all higher 
order terms results in the following set of linear equations. 


A P 2 (fc) ' 


8P 2 ( fc ) 
dS 2 

dP 2 ( fc ) 
dS n 

8P ?I ( fc ) 
d\V 2 \ 

8P ? , ( fc ) 
d\V n \ 


2c 

<1 

A Pfi k) 


dPn ( fc ) 

as 2 

8P n ( k ) 
' dS n 

8P n ( fc ) 

mi 

8P n ( fc ) 
d\Vn\ 


A 5 { n k) 

O* •' 
<1 


dQ 2 ( fc ) 
8S 2 

8Q 2 ( fc ) 

' 8S n 

8Q 2 ( fe ) 
d\V 2 \ 

aq 2 ( fc ) 

8\V„\ 


A|F 2 (fc) | 

A Q ( n k) _ 


8Q n W 
då 2 

8Qn W 

■ dSn 

8Qn W 

8m 

dQn ( fc ) 

»|V„| J 


A|Fi fc) | 


In the above equation, bus 1 is assumed to be the slack bus. The Jacobian matrix 
gives the linearized relationship between small changes in voltage angle Ad) ; 
and voltage magnitude | with the small changes in real and reactive power 

A and A Q\ k \ Elements of the Jacobian matrix are the partial derivatives of 

(6.52) and (6.53), evaluated at A and A|V^|. In short form, it can be written 
as 


A P 


Ji 

Ji 

' AS ' 

AQ 


. Js 

J 4 

. A|V| . 


For voltage-controlled buses, the voltage magnitudes are known. Therefore, if m 
buses of the system are voltage-controlled, m equations involving AQ and AF 
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and the corresponding columns of the Jacobian matrix are eliminated. Accordingly, 
there are n — 1 real power constraints and n — 1 — m reactive power constraints, 
and the Jacobian matrix is of order (2 n — 2 — m) x (2 n — 2 — m). J\ is of the 
order (n — 1 ) x (n — 1 ), J 2 is of the order (n — 1 ) x (n — 1 — m), J 3 is of the 
order (n — 1 — m) x (n — 1), and J 4 is of the order (n — 1 — m) x (n — 1 — m). 

The diagonal and the off-diagonal elements of J\ are 


= ^Z\V i \\V j \\Y ij \sm(e ij -Si + 5 j ) 

= -MIW« I s in(% -4 + Sj) j yé i 


(6.55) 


(6.56) 


The diagonal and the off-diagonal elements of J 2 are 


2\V i \\Y ii \ cos ØH + IVjWYijl cosiøij — Si + 6j) (6.57) 


Q|y.| — | Vi 1 1 | COS (Øj j Si + Sj) j^i 

The diagonal and the off-diagonal elements of J 3 are 


(6.58) 


= E \ViWVjl\Yij\ cos (% - Si + Sj) (6.59) 

1 

= -1^11^11^1 cos(% -5i + Sj) j yé i (6.60) 

The diagonal and the off-diagonal elements of J 4 are 


-2|Vi||y«| sin 9a - £ \Vj\\Y i5 1 sin(% - % + Sj) (6.61) 


= -\Vi\\Yij | sm(9ij ~Si + Sj) j / i 


(6.62) 


The terms A and A are the difference between the scheduled and calcu- 
lated values, known as the power residuals, given by 

AP/ fc) = Pf ch - jf } (6.63) 

AQ\ k) = Q\ ch - Q {k) (6.64) 
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The new estimates for bus voltages are 

+1) = S\ k) + Atff ] (6.65) 

|v/ fc+1) | = |F/ fc) | + A|V; (fc) | (6.66) 

The procedure for power flow solution by the Newton-Raphson method is as 
follows: 

1. For load buses, where P? ch and Qf ch are specified, voltage magnitudes and 

phase angles are set equal to the slack bus values, or 1.0 and 0.0, i.e., | = 

1.0 and sj°^ — 0.0. For voltage-regulated buses, where |Vj| and Pf ch are 
specified, phase angles are set equal to the slack bus angle, or 0, i.e., = 0. 

2. For load buses, and are calculated from (6.52) and (6.53) and 
AP{' k> and A are calculated from (6.63) and (6.64). 

3. For voltage-controlled buses, P-^ and A are calculated from (6.52) and 
(6.63), respectively. 

4. The elements of the Jacobian matrix (J\, J 2 , J3, and J4) are calculated 
from (6.55) - (6.62). 

5. The linear simultaneous equation (6.54) is solved directly by optimally or- 
dered triangular factorization and Gaussian elimination. 

6. The new voltage magnitudes and phase angles are computed from (6.65) and 

( 6 . 66 ). 

7. The process is continued until the residuals A P^ and A are less than 
the specified accuracy, i.e., 

|AP/ fc) | < e 

|AQf } |<e (6.67) 

The power flow solution by the Newton-Raphson method is demonstrated in 
the following example. 

Example 6.10 

Obtain the power flow solution by the Newton-Raphson method for the system of 
Example 6.8. 
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Line impedances converted to admittances are j/i 2 = 10 — j20, 2/13 = 10 — .7 30, 
and 1/23 = 16 — j32. This results in the bus admittance matrix 


Ybus — 


20 - j50 
— 10 + y20 
-10 + y30 


-10 + j20 
26 - j52 
-16 + j32 


-10 + j30 * 
-16 + j32 
26 - j62 


Converting the bus admittance matrix to polar form with angles in radian yields 


Ybus — 


53.85165Z -1.9029 
22.36068Z2.0344 
31. 62278Z 1.8925 


22.36068 Z 2.0344 
58.13777Z-1.1071 
35.77709Z2.0344 


31.62278Z 1.8925 ' 
35.77709Z 2.0344 
67.23095Z -1.1737 


From (6.52) and (6.53), the expressions for real power at bus 2 and 3 and the 
reactive power at bus 2 are 

P 2 = | ^ 2 1 1 ^ 1 1 1^21 1 COs (#21 — <*>2 + h) + | V 2 1 1^22 1 COS $22 + 

IV2IIV3KY23I COs($ 23 — 82 + 83) 

P 3 = |V 3 ||Vi||ni|cOB (031 - 5 3 + * 1 ) + l^ll^linal cos (032 - 

83 + 82) + S V3 2 1 1+33 1 COS 633 

Q 2 = — |IZ 2 ||Fi||Y 2 i|sin(02i - J 2 + J a ) - |V 2 2 ||Y 22 |sin022 - 
!^2||Y 3 ||Y 2 3|si n (023- J2 + J3) 

Elements of the Jacobian matrix are obtained by taking partial derivatives of the 
above equations with respect to J 2 , J 3 and |Y 2 |. 

1^21 1 Ya ||Y 2 i I sin(0 21 - å 2 + 8 X ) + I^IIVall^l 
sin (#23 - 8 2 + 83) 

-|V r J||V 3 ||y 23 |sin(0 2 3-^ + « 3 ) 

|Fi||^2i!cos(e 2 i - J 2 + Ji) + 2|F 2 ||Y 22 |cos0 22 + 

|^||^ 23 | COS ( 0 23 -82 + 83) 

— I^ll^ll^bl sin(032 — 83 + 82) 

|V 3 ||Vi||y 3 i|sm(fl 31 - J 3 + h) + |V3||y 2 ||i32| 

sin (^ 32 — 83 + J 2 ) 

|Y 3 (|Y 3 2| COs (^32 — J3 + 82) 


dP 2 

dS 2 

dP 2 

dS z 

dP 2 

d\V 2 \ 

dP z 

dS 2 

dP z 

d8 z 

dP z 
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~ = I^M^||y 2 i|cos(ø 21 - 6 2 + s 1 ) + \V 2 m\YM\ 

COS (023 - $2 + 5s) 

—■ = -|^||V3i|y 23 |c 0 s(023 - ^2 + Sz) 

= — | V]. 1 1^21 1 sin (0 2 i - <5 2 + S x ) - 2|y 2 ||y 2 2|sin0 2 2 - 
1^31 |V23 1 sin (023 - $2 + S 3 ) 

The load and generation expressed in per units are 

csch _ (400 + j‘250) , n , 0 _ 


gsch. _ 


= —4.0 — j2.5 pu 


Pi* = ^ = 2.° P” 

The slack bus voltage is Vi = 1.05Z0 pu, and the bus 3 voltage magnitude is 
j V 3 1 = 1.04 pu. Starting with an initial estimate of |V 2 0 ^| = 1.0, 4°^ = 0.0, and 
63 ® = 0.0, the power residuals are computed from (6.63) and (6.64) 

A P 2 (0) = Pi ch - P 2 (0) = -4.0 - (-1.14) - -2.8600 
A P 3 (0) = P| cft - P 3 (0) = 2.0 - (0.5616) = 1.4384 


A Qf 


psch _ pw _ 2.0 - (0.5616) = 1.4384 
Q 2 h ~ Q? = -2-5 - (-2.28) = - 0.2200 


Evaluating the elements of the Jacobian matrix with the initial estimate, the set of 
linear equations in the first iteration becomes 


-2.8600 

1.4384 

- 0.2200 


54.28000 -33.28000 24.86000 

-33.28000 66.04000 -16.64000 

-27.14000 16.64000 49.72000 


A|V 2 (0) | 


Obtaining the solution of the above matrix equation, the new bus voltages in the 
first iteration are 


A4 0) = -0.045263 4 1} = 0 + (-0.045263) = -0.045263 

A4 0) = -0.007718 4 1} = 0 + (-0.007718) = -0.007718 

A|V 2 (0) | =-0.026548 jV 2 {1) | = 1 + (-0.026548) = 0.97345 

Voltage phase angles are in radians. For the second iteration, we have 


-0.099218 


51.724675 -31.765618 21.302567 


0.021715 = -32.981642 65.656383 -15.379086 


a4 1} 


-0.050914 -28.538577 17.402838 48.103589 I I A|V 2 (1) | 
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and 

a4 1} = -0.001795 4 2) = -0.045263 + (-0.001795) = -0.04706 

a4 1} = -0.000985 4 2) = -0.007718 + (-0.000985) = -0.00870 

A|V 2 (1) | = -0.001767 |V 2 (2) | = 0.973451 + (-0.001767) = 0.971684 


For the third iteration, we have 


' -0.000216 ‘ 


51.596701 -31.693866 21.147447 ' 


a4 2) 

0.000038 

= 

-32.933865 65.597585 -15.351628 


a4 2) 

-0.000143 


-28.548205 17.396932 47.954870 


. A|V 2 (2) | _ 


and 

A 4 2) = -0.000038 4 3) = -0.047058 + (-0.0000038) = -0.04706 

A4 2) = -0.0000024 4 3) = -0.008703 + (-0.0000024) = 0.008705 

A|V^ 2) j = -0.0000044 |V 2 (3) | = 0.971684 + (-0.0000044) = 0.97168 

The solution converges in 3 iterations with a maximum power mismatch of 2.5 x 
10~ 4 with V 2 = 0.97168Z -2.696° and V 3 = 1.04Z-0.4988 0 . From (6.52) and 
(6.53), the expressions for reactive power at bus 3 and the slack bus real and reac- 
tive powers are 

Qz = — iV 3 ||Vi||l^i| sin(6> 3 i - <5 3 + S x ) - IF 3 IIW 32 I 

sin(032 — S 3 + J 2 ) — | Vz ] 2 II 33 1 sin #33 

Pi = |^i| 2 |ni|cosØn + |^ 1 ||^ 2 ||yi 2 |cos( 6» 1 2 - <5i + <5 2 ) + | Vi 1 1 V 3 1 

IV 13 I COS (013 — ^1 + 53 ) 

Qi = -IVxplFnlsinØn - |Vi||V a ||yi2|sin(0 12 - S x + <5 2 ) - |V X ||V 3 | 
|Vi 3 |sin(Øi 3 - 5i +<J 3 ) 

Upon substitution, we have 

Q 3 = 1.4617 pu 
Pi = 2.1842 pu 
Qi = 1.4085 pu 

Finally, the line flows are calculated in the same manner as the line flow calcula- 
tions in the Gauss-Seidel method described in Example 6.7, and the power flow 
diagram is as shown in Figure 6.13. 

A program named Ifnewton is developed for power flow solution by the 
Newton-Raphson method for practical power systems. This program must be pre- 
ceded by the Ifybus program, busout and lineflow programs can be used to print 
the load flow solution and the line flow results. The format is the same as the 
Gauss-Seidel. The following is a brief description of the Ifnewton program. 
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Ifnewton This program obtains the power flow solution by the Newton-Raphson 
method and requires the busdata and the linedata files described in Sec- 
tion 6.9. It is designed for the direct use of load and generation in MW and 
Mvar, bus voltages in per unit, and angle in degrees. Loads and generation 
are converted to per unit quantities on the base MVA selected. A provision 
is made to maintain the generator reactive power of the voltage-controlled 
buses within their specified limits. The violation of reactive power limit may 
occur if the specified voltage is either too high or too low. In the second it- 
eration, the var calculated at the generator buses are examined. If a limit is 
reached, the voltage magnitude is adjusted in steps of 0.5 percent up to ±5 
percent to bring the var demand within the specified limits. 


Example 6.11 

Obtain the power flow solution for the IEEE-30 bus test system by the Newton- 
Raphson method. 

The data required is the same as in Example 6.9 with the following commands 

clear 7. clears all variables from the workspace. 

basemva = 100; accuracy = 0.001; maxiter = 12; 

busdata = [ same as in Example 6.9]; 

linedata = [ same as in Example 6.9]; 

lfybus 7, Forms the bus admittance matrix 

Ifnewton 7. Power flow solution by Newton-Raphson method 
busout 7« Prints the power flow solution on the screen 

lineflow 7. Computes and displays the line flow and losses 


The output of Ifnewton is 


Power Flow Solution by Newton-Raphson Method 
Maximum Power mismatch = 7.54898e-07 
No. of iterations = 4 


Bus Voltage Angle Load — Generation — Injected 


No . 

Mag. 

Degree 

MW 

Mvar 

MW 

Mvar 

Mvar 

1 

1.060 

0.000 

0.000 

0.000 

260.998 

-17.021 

0.00 

2 

1.043 

-5.497 

21.700 

12.700 

40.000 

48.822 

0.00 

3 

1.022 

-8.004 

2.400 

1.200 

0.000 

0.000 

0.00 
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4 

1 . 

013 

-9. 

661 

7. 

600 

1 . 

600 

0 . 

000 

0 . 

000 

0 . 

00 

5 

1 . 

010 

-14. 

,381 

94. 

200 

19. 

000 

0 . 

000 

35. 

975 

0 . 

00 

6 

1. 

,012 

-11. 

398 

0 . 

000 

0. 

000 

0 . 

000 

0 . 

000 

0 . 

00 

7 

1 . 

.003 

-13. 

,150 

22. 

800 

10. 

900 

0 . 

000 

0 . 

000 

0 . 

00 

8 

1 . 

.010 

-12. 

.115 

30. 

000 

30. 

000 

0 . 

000 

30. 

826 

0 . 

00 

9 

1 . 

.051 

-14. 

,434 

0 . 

000 

0 . 

000 

0 . 

000 

0 . 

000 

0 . 

00 

10 

1 . 

.044 

-16. 

,024 

5. 

800 

2. 

000 

0 . 

000 

0 . 

000 

19. 

,00 

11 

1 . 

.082 

-14. 

.434 

0 . 

000 

0 . 

000 

0 . 

000 

16. 

119 

0 . 

.00 

12 

1 . 

.057 

-15. 

,302 

11. 

200 

7. 

,500 

0 . 

000 

0 . 

.000 

0 . 

.00 

13 

1 . 

.071 

-15. 

,302 

0 . 

000 

0. 

,000 

0 . 

,000 

10. 

,423 

0 . 

.00 

14 

1 . 

.042 

-16. 

.191 

6. 

200 

1 . 

.600 

0 . 

000 

0 . 

,000 

0. 

.00 

15 

1 . 

.038 

-16. 

.278 

8. 

200 

2. 

.500 

0 . 

,000 

0 . 

,000 

0 . 

.00 

16 

1 . 

.045 

-15. 

.880 

3. 

500 

1 . 

.800 

0 . 

.000 

0 . 

.000 

0 . 

.00 

17 

1 . 

.039 

-16. 

.188 

9. 

000 

5, 

.800 

0. 

,000 

0 . 

.000 

0. 

.00 

18 

1 , 

.028 

-16. 

.884 

3. 

,200 

0. 

.900 

0. 

.000 

0 . 

.000 

0 . 

.00 

19 

1 . 

.025 

-17. 

.052 

9. 

500 

3. 

.400 

0 . 

.000 

0 , 

.000 

0. 

.00 

20 

1 . 

.029 

-16. 

.852 

2. 

.200 

0 . 

.700 

0 

.000 

0, 

.000 

0 . 

.00 

21 

1 

.032 

-16. 

.468 

17. 

.500 

11, 

.200 

0 , 

.000 

0 . 

.000 

0 

.00 

22 

1, 

.033 

-16, 

.455 

0, 

.000 

0 

.000 

0 

.000 

0, 

.000 

0, 

.00 

23 

1 

.027 

-16 

.662 

3. 

.200 

1, 

.600 

0 

.000 

0 

.000 

0 

.00 

24 

1 

.022 

-16 

.830 

8, 

.700 

6 

.700 

0 

.000 

0 

.000 

4 

.30 

25 

1 

.019 

-16 

.424 

0 

.000 

0 

.000 

0 

.000 

0 

.000 

0 

.00 

26 

1 

.001 

-16 

.842 

3, 

.500 

2 

.300 

0 

.000 

0 

.000 

0 

.00 

27 

1 

.026 

-15 

.912 

0 

.000 

0 

.000 

0 

.000 

0 

.000 

0 

.00 

28 

1 

.011 

-12 

.057 

0 

.000 

0 

.000 

0 

.000 

0 

.000 

0 

.00 

29 

1 

.006 

-17 

.136 

2 

.400 

0 

.900 

0 

.000 

0 

.000 

0 

.00 

30 

0 

.995 

-18 

.015 

10 

.600 

1 

.900 

0 

.000 

0 

.000 

0 

.00 

Total 




283 

.400 

126 

.200 

300 

.998 

125 

.144 

23 

.30 


The output of the lineflow is the same as the line flow output of Example 6.9 with 
the power mismatch as dictated by the Newton-Raphson method. 


6.11 FAST DECOUPLED POWER FLOW SOLUTION 

Power system transmission lines have a very high X/R ratio. For such a system, 
real power changes A P are less sensitive to changes in the voltage magnitude 
and are most sensitive to changes in phase angle Ad. Similarly, reactive power is 
less sensitive to changes in angle and are mainly dependent on changes in voltage 
magnitude. Therefore, it is reasonable to set elements J 2 and J 3 of the Jacobian 
matrix to zero. Thus, (6.54) becomes 


Jl 

0 

' A<5 ■ 

L 0 

J 4 . 

. A|V| . 


AP 

A Q 


( 6 . 68 ) 
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dP 

AP = J 1 AS = [— -] AS 
od 

AQ = J a A |V| = Æ] A|V| 


(6.69) 


(6.70) 


(6.69) and (6.70) show that the matrix equation is separated into two decoupled 
equations requiring considerably less time to solve compared to the time required 
for the solution of (6.54). Furthermore, considerable simplification can be made to 
eliminate the need for recomputing J\ and J 4 during each iteration. This procedure 
results in the decoupled power flow equations developed by Stott and Alsac[75- 
76]. The diagonal elements of J\ described by (6.55) may be written as 

BP n 

7^ = E IVill^ll^lsin^ - Si + Sj) - \Vi\ 2 \Y«\éme u 

OOi j=1 

Replacing the first term of the above equation with —Qi, as given by (6.53), results 


BP- 

^ = -Q i -N 2 |y ii |sinø ii 
= -Qi-\V i \ 2 B ii 

Where Bu = |Yj;| sin du is the imaginary part of the diagonal elements of the bus 
admittance matrix. Bu is the sum of susceptances of all the elements incident to bus 
i. In a typical power system, the self-susceptance Bu Qi, and we may neglect 
Qi. Further simplification is obtained by assuming ]Vj | 2 æ |V*|, which yields 


= -\Vi\Bu 


(6.71) 


Under normal operating conditions, 5j — di is quite small. Thus, in (6.56) assuming 
Ou — 5i + 6j æ Ou, the off-diagonal elements of J\ becomes 

BP 

= -\ViWVi\By 


Further simplification is obtained by assuming \Vj\ « 1 

S = -IKI B l} 
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Similarly, the diagonal elements of J 4 described by (6.61) may be written as 

BO- n 

^j~| = -M|y«| Sin du - £ miVjWYijl sin(% - * + Sj) 

replacing the second term of the above equation with -Qi, as given by (6.53), 
results in 

= — IViUYnl sinØji + Qi 

Again, since Bu = Yu sin du > Q it Qi may be neglected and (6.61) reduces to 

m = - mBu (6 - 73) 

Likewise in (6.62), assuming Oij — Si + 6j ~ Oij yields 

m = (6 - 74) 

With these assumptions, equations (6.69) and (6.70) take the following form 

A P 

7777 = -B' M (6.15) 


= —B" A|V| 


(6.76) 


Here, B 1 and B" are the imaginary part of the bus admittance matrix Y& us . Since 
the elements of this matrix are constant, they need to be triangularized and in- 
verted only once at the beginning of the iteration. B' is of order of (n — 1). For 
voltage-controlled buses where | 1 ^| and Pi are specified and Qi is not specified, 
the corresponding row and column of Yb ua are eliminated. Thus, B" is of order of 
(n — 1 — m), where m is the number of voltage-regulated buses. Therefore, in the 
fast decoupled power flow algorithm, the successive voltage magnitude and phase 
angle changes are 



(6.77) 


(6.78) 


The fast decoupled power flow solution requires more iterations than the Newton- 
Raphson method, but requires considerably less time per iteration, and a power 
flow solution is obtained very rapidly. This technique is very useful in contingency 
analysis where numerous outages are to be simulated or a power flow solution is 
required for on-line Control. 
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Example 6.12 

Obtain the power flow solution by the fast decoupled method for the system of Ex- 
ample 6 . 8 . 


The bus admittance matrix of the system as obtained in Example 6.10 is 


ms — 


20 — y'50 -lO + j/20 -10 + j30' 
-10 + j20 26 — j52 — 16 + j32 

-10 + j30 -16 + j32 26 - j 62 . 


In this system, bus 1 is the slack bus and the corresponding bus susceptance matrix 
for evaluation of phase angles Afø and A(5 3 is 


B' = 



32 ' 
-62 


The inverse of the above matrix is 


[BT 1 


-0.028182 -0.014545 ' 
-0.014545 -0.023636 


From (6.52) and (6.53), the expressions for real power at bus 2 and 3 and the 
reactive power at bus 2 are 

P 2 = |Vb||Vi||y 2 i|cos(ø 21 - <5 2 + Si) + |y 2 2 ||r 22 |cosø 22 
+ |^||V3||y23|00B(023-5 2 + 5s) 

P 3 = |V 3 ||Vi||Ysi| cos(^ 3 i — <^3 + <5i) + I V 3 1| V 2 1 IV 32 1 cos(é? 3 2 
-£3 + * 2 ) + |V 3 2 ||y 3 3lc 0 S033 
Q 2 = — | V 2 1 1 Fi 1 1 V 21 i sin(^ 2 i. - <52 + Si) - |V 2 2 ||y 2 2|sm022 
— I ^ 2 1 1 1 1 ^23 1 sin(^23 — <52 + £ 3 ) 


The load and generation expressed in per units are 


osch 

*- > 2 — 


P i sch 

3 — 


(400 + j'250) 


200 

100 


100 

= 2.0 pu 


-4.0 - j 2.5 pu 


The slack bus voltage is Vi = 1.05Z0 pu, and the bus 3 voltage magnitude is 
| V 3 I = 1.04 pu. Starting with an initial estimate of |V 2 °^| = 1.0, <5 2 °^ = 0.0, and 
4 0) = 0.0, the power residuals are computed from (6.63) and (6.64) 


AP 2 (0) = Pi ch - P 2 (0) = -4.0 - (-1.14) = -2.86 
A P 3 (0) = Pi ch - P 3 (0) = 2.0 - (0.5616) = 1.4384 
AQ ( 2 0) = Q s 2 ch - q¥ ] = -2.5 - (-2.28) = -0.22 
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The fast decoupled power flow algorithm given by (6.77) becomes 


‘ a4°> ' 

• 

A \ 



-0.028182 -0.014545 
-0.014545 -0.023636 



f -2.8600 1 



|£m 

iwb 

= 

• 

L 1-04 J 

- 


-0.060483 

-0.008909 


Since bus 3 is a regulated bus, the corresponding row and column of B' are elimi- 
nated and we get 


From (6.78), we have 


B" = [-52] 


A|F 2 | 


- 1 ' 


— . 22 ' 

.52. 


. 10 . 


-0.0042308 


The new bus voltages in the first iteration are 


a4 0) = -0.060483 
A4 0) = -0.008989 
A|F 2 (0) | = -0.0042308 


4 1} = 0 + (-0.060483) = -0.060483 
4 1} = 0 + (-0.008989) = -0.008989 
|y 2 (1) | = 1 + (-0.0042308) = 0.995769 


The voltage phase angles are in radians. The process is continued until power resid- 
uals are within a specified accuracy. The result is tabulated in the table below. 


Iter 

$2 

<>3 

m 

A P 2 

ap 3 

A Q 2 

1 

-0.060482 

-0.008909 

0.995769 

-2.860000 

1.438400 

- 0.220000 

2 

-0.056496 

-0.007952 

0.965274 

0.175895 

-0.070951 

-1.579042 

3 

-0.044194 

-0.008690 

0.965711 

0.640309 

-0.457039 

0.021948 

4 

-0.044802 

-0.008986 

0.972985 

-0.021395 

0.001195 

0.365249 

5 

-0.047665 

-0.008713 

0.973116 

-0.153368 

0.112899 

0.006657 

6 

-0.047614 

-0.008645 

0.971414 

0.000520 

0.002610 

-0.086136 

7 

-0.046936 

-0.008702 

0.971333 

0.035980 

-0.026190 

-0.004067 

8 

-0.046928 

-0.008720 

0.971732 

0.000948 

-0.001411 

0.020119 

9 

-0.047087 

-0.008707 

0.971762 

-0.008442 

0.006133 

0.001558 

10 

-0.047094 

-0.008702 

0.971669 

-0.000470 

0.000510 

-0.004688 

11 

-0.047057 

-0.008705 

0.971660 

0.001971 

-0.001427 

-0.000500 

12 

-0.047054 

-0.008706 

0.971681 

0.000170 

-0.000163 

0.001087 

13 

-0.047063 

-0.008706 

0.971684 

-0.000458 

0.000330 

0.000151 

14 

-0.047064 

-0.008706 

0.971680 

-0.000053 

0.000048 

-0.000250 


Converting phase angles to degrees the final solution is V 2 = 0.97168Z— 2.696° 
and V 3 = 1.04Z— 0.4988°. Using (6.52) and (6.53) as in Example 6.10, the reactive 
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power at bus 3 and the slack bus real and reactive powers are 

Q 3 = 1.4617 pu 
Pi = 2.1842 pu 
Q\ = 1.4085 pu 

The fast decoupled power flow for this example has taken 14 iterations with the 
maximum power mismatch of 2.5 x 10~ 4 pu compared to the Newton-Raphson 
method which took only three iterations. The highest X/R ratio of the transmission 
lines in this example is 3. For systems with a higher X/R ratio, the fast decoupled 
power flow method converges in relatively fewer iterations. However, the number 
of iterations is a function of system size. 

Finally, the line flows are calculated in the same manner as the line flow cal- 
culations in the Gauss-Seidel method described in Example 6.7, and the power flow 
diagram is as shown in Figure 6.13. 

A program named decouple is developed for power flow solution by the fast 
decoupled method for practical power systems. This program must be preceded by 
the Ifybus program, busout and lineflow programs can be used to print the load 
flow solution and the line flow results. The format is the same as the Gauss-Seidel 
method. The following is a brief description of the decouple program: 


decouple This program finds the power flow solution by the fast decouple method 
and requires the busdata and the linedata files described in Section 6.9. It is 
designed for the direct use of load and generation in MW and Mvar, bus volt- 
ages in per unit, and angle in degrees. Loads and generation are converted 
to per unit quantities on the base MVA selected. A provision is made to 
maintain the generator reactive power of the voltage-controlled buses within 
their specified limits. The violation of reactive power limit may occur if the 
specified voltage is either too high or too low. In the lOth iteration, the vårs 
calculated at the generator buses are examined. If a limit is reached, the volt- 
age magnitude is adjusted in steps of 0.5 percent up to ±5 percent to bring 
the var demand within the specified limits. 


Example 6.13 

Obtain the power flow solution for the IEEE-30 bus test system by the fast decou- 
pled method. 


Data required is the same as in Example 6.9 with the following commands 


246 6. POWER FLOW ANALYSIS 


clear 7, clears all variables from the vorkspace. 

basemva = 100; accuracy = 0.001; maxiter = 20; 

busdata = [ same as in Example 6.9 ]; 

linedata = [ same as in Example 6.9]; 

lfybus % Forms the bus admittance matrix 

decouple 7, Power flow solution by fast decoupled method 

busout I, Prints the power flow solution on the screen 

lineflow 7« Computes and displays the line flow and losses 

The output of decouple is 

Power Flow Solution by Fast Decoupled Method 
Maximum Power mismatch = 0.000919582 
No. of iterations = 15 


Bus Voltage Angle Load — Generation — Injected 


No . 

Mag. 

Degree 

MW 

1 

1 

.060 

0 . 

.000 

0 , 

.000 

2 

1 

.043 

-5, 

.497 

21. 

.700 

3 

1 

.022 

-8. 

.004 

2, 

.400 

4 

1 

.013 

-9. 

.662 

7. 

.600 

5 

1 

.010 

-14, 

.381 

94, 

.200 

6 

1 

.012 

-11. 

,398 

0. 

.000 

7 

1 

.003 

-13, 

.149 

22, 

.800 

8 

1 

.010 

-12, 

.115 

30, 

.000 

9 

1 

.051 

-14. 

,434 

0 , 

.000 

10 

1 

.044 

-16. 

.024 

5. 

.800 

11 

1 

.082 

-14. 

.434 

0. 

.000 

12 

1 

.057 

-15. 

,303 

11, 

.200 

13 

1, 

.071 

-15. 

.303 

0 , 

.000 

14 

1 

.042 

-16. 

,198 

6, 

.200 

15 

1 

.038 

-16. 

,276 

8. 

.200 

16 

1 . 

.045 

-15. 

,881 

3, 

.500 

17 

1 

.039 

-16. 

.188 

9. 

.000 

18 

1 . 

.028 

-16. 

,882 

3, 

.200 

19 

1. 

.025 

-17. 

,051 

9. 

.500 

20 

1 . 

.029 

-16. 

.852 

2. 

.200 

21 

1. 

.032 

-16. 

,468 

17. 

.500 

22 

1. 

.033 

-16. 

,454 

0 . 

.000 

23 

1 . 

.027 

-16. 

.661 

3. 

.200 

24 

1 , 

.022 

-16. 

,829 

8, 

.700 

25 

1 . 

.019 

-16. 

,423 

0. 

,000 

26 

1 . 

.001 

-16. 

840 

3. 

.500 


Mvar MW Mvar Mvar 


0 , 

.000 

260 

.998 

-17, 

.021 

0 . 

.00 

12 

.700 

40, 

.000 

48, 

.822 

0 , 

.00 

1 

.200 

0, 

.000 

0 , 

.000 

0 . 

.00 

1 , 

.600 

0 , 

.000 

0 

,000 

0 , 

,00 

19. 

.000 

0 , 

.000 

35, 

.975 

0 , 

.00 

0 . 

.000 

0 , 

.000 

0 . 

.000 

0 . 

.00 

10. 

.900 

0 . 

.000 

0 , 

.000 

0 . 

.00 

30. 

.000 

0 , 

.000 

30, 

.828 

0. 

.00 

0 . 

.000 

0 , 

.000 

0 , 

.000 

0 . 

.00 

2. 

.000 

0 , 

.000 

0 

.000 

19. 

.00 

0 . 

.000 

0 , 

.000 

16. 

.120 

0 . 

.00 

7. 

.500 

0 , 

.000 

0 

.000 

0 . 

.00 

0 . 

.000 

0 . 

.000 

10. 

.421 

0. 

.00 

1 . 

.600 

0, 

.000 

0 , 

.000 

0 . 

.00 

2. 

.500 

0 , 

.000 

0 , 

.000 

0 , 

.00 

1 , 

.800 

0 . 

,000 

0 , 

.000 

0 , 

.00 

5, 

.800 

0 . 

.000 

0 

.000 

0 , 

.00 

0 . 

.900 

0. 

.000 

0 , 

.000 

0 . 

.00 

3, 

.400 

0. 

.000 

0 , 

.000 

0 . 

.00 

0 . 

.700 

0 . 

.000 

0 , 

.000 

0 , 

.00 

11, 

.200 

0. 

,000 

0 , 

.000 

0, 

.00 

0 , 

.000 

0 , 

.000 

0 

.000 

0 , 

.00 

1, 

.600 

0 . 

.000 

0 , 

.000 

0. 

.00 

6. 

.700 

0 , 

.000 

0 , 

.000 

4, 

.30 

0. 

.000 

0. 

,000 

0 , 

.000 

0 , 

.00 

2. 

.300 

0 . 

.000 

0, 

.000 

0, 

.00 
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27 1 . 026 

- 15.912 

0.000 

0.000 

0.000 

0.000 

0.00 

28 1.011 

- 12.057 

0.000 

0.000 

0.000 

0.000 

0.00 

29 1 . 006 

- 17.136 

2.400 

0.900 

0.000 

0.000 

0.00 

30 0.995 

- 18.014 

10.600 

1.900 

0.000 

0.000 

0.00 

Total 


283.400 

126.200 

300.998 

125.145 

23.30 


The output of the lineflow is the same as the line flow output of Example 6.9 with 
the power mismatch as dictated by the fast decoupled method. 

PROBLEMS 

6.1. A power system network is shown in Figure 6.17. The generators at buses 
1 and 2 are represented by their equivalent current sources with their reac- 
tances in per unit on a 100-MVA base. The lines are represented by 7r model 
where series reactances and shunt reactances are also expressed in per unit 
on a 100 MYA base. The loads at buses 3 and 4 are expressed in MW and 
Mvar. 

(a) Assuming a voltage magnitude of 1.0 per unit at buses 3 and 4, convert 
the loads to per unit impedances. Convert network impedances to admit- 
tances and obtain the bus admittance matrix by inspection. 

(b) Use the function Y = ybus(zdata) to obtain the bus admittance matrix. 
The function argument zdata is a matrix containing the line bus numbers, 
resistance and reactance. (See Example 6.1.) 



FIGURE 6.17 

One-line diagram for Problem 6.1. 
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6.2. A power system network is shown in Figure 6.18. The values marked are 
impedances in per unit on a base of 100 MVA. The currents entering buses 1 
and 2 are 

h = 1.38 — y'2.72 pu 
I 2 = 0.69 — jl.36 pu 

(a) Determine the bus admittance matrix by inspection. 

(b) Use the function Y = ybus(zdata) to obtain the bus admittance matrix. 
The function argument zdata is a matrix containing the line bus numbers, 
resistance and reactance. (See Example 6.1.) Write the necessary MATLAB 
commands to obtain the bus voltages. 


0 



6 


FIGURE 6.18 

One-line diagram for Problem 6.2. 

6.3. Use Gauss-Seidel method to find the solution of the following equations 

x\ + X 1 X 2 = 10 
Xi+ X 2 = 6 

with the following initial estimates 

(a) xf^ = 1 and = 1 

(b) = 1 and x^ = 2 

Continue the iterations until |Axj^| and lAx^l are less than 0.001. 
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6.4. A fourth-order polynomial equation is given by 

x 4 - 21x 3 -i- U7x 2 - 379x + 252 = 0 

(a) Use Newton-Raphson method and hånd calculations to find one of the 
roots of the polynomial equation. Start with the initial estimate of x = 0 
and continue until |Aa4 fc )| < 0.001. 

(b) Write a MATLAB program to find the roots of the above polynomial by 
Newton-Raphson method. The program should prompt the user to input the 
initial estimate. Run using the initial estimates of 0, 3, 6, 10. 

(c) Check your answers using the MATLAB function r = roots(A), where A 
is a row vector containing the polynomial coefficients in descending powers. 

6.5. Use Newton-Raphson method and hånd calculation to find the solution of 
the following equations: 

x^ — 2xi — X 2 = 3 
x\ + = 41 

(a) Start with the initial estimates of = 2, x^ — 3 . Perform three itera- 
tions. 

(b) Write a MATLAB program to find one of the solutions of the above equa- 
tions by Newton-Raphson method. The program should prompt the user to 
input the initial estimates. Run the program with the above initial estimates. 

6.6. In the power system network shown in Figure 6.19, bus 1 is a slack bus with 
Vi = 1.0Z0° per unit and bus 2 is a load bus with S 2 = 280 MW + j60 
Mvar. The line impedance on a base of 100 MVA is Z — 0.02 + j0.04 per 
unit. 

(a) Using Gauss-Seidel method, determine V<i ■ Use an initial estimate of 

= 1.0 + j 0.0 and perform four iterations. 

(b) If after several iterations voltage at bus 2 converges to V 2 = 0.90— jO. 10, 
determine S\ and the real and reactive power loss in the line. 



FIGURE 6.19 

One-line diagram for Problem 6.6. 
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FIGURE 6.20 

One-line diagram for Problem 6.7. 


6.7. Figure 6.20 shows the one-line diagram of a simple three-bus power system 
with generation at bus 1. The voltage at bus 1 is Vi = 1.0/0° per unit. 
The scheduled loads on buses 2 and 3 are marked on the diagram. Line 
impedances are marked in per unit on a 100-MVA base. For the purpose 
of hånd calculations, line resistances and line charging susceptances are ne- 
glected. 

(a) Using Gauss-Seidel method and initial estimates of = 1.0 + jO and 

V 3 ^ = 1.0 + jO, determine V2 and V3. Perform two iterations. 

(b) If after several iterations the bus voltages converge to 

V2 = 0.90 — y'0.10 pu 
V 3 = 0.95 - j0.05 pu 

determine the line flows and line losses and the slack bus real and reactive 
power. Construct a power flow diagram and show the direction of the line 
flows. 

(c) Check the power flow solution using the Ifgauss and other required pro- 
grams. (Refer to Example 6.9.) Use a power accuracy of 0.00001 and an 
acceleration factor of 1.0. 

6 .8. Figure 6.21 shows the one-line diagram of a simple three-bus power system 
with generation at buses 1 and 3. The voltage at bus 1 is Vi = 1.025/0° per 
unit. Voltage magnitude at bus 3 is fixed at 1.03 pu with a real power gener- 
ation of 300 MW. A load consisting of 400 MW and 200 Mvar is taken from 
bus 2. Line impedances are marked in per unit on a 100-MVA base. For the 
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FIGURE 6.21 

One-line diagram for Problem 6.8. 

purpose of hånd calculations, line resistances and line charging susceptances 
are neglected. 

(a) Using Gauss-Seidel method and initial estimates of = 1.0 + j0 and 

V 3 ^ = 1.03 + jO and keeping j V 3 I = 1.03 pu, determine the phasor values 
of V 2 and V 3 . Perform two iterations. 

(b) If after several iterations the bus voltages converge to 

V 2 = 1.001243Z— 2.1° = 1.000571 - jO.0366898 pu 
y 3 = 1.03Z 1.36851° = 1.029706 + jO.0246 pu 

determine the line flows and line losses and the slack bus real and reactive 
power. Construct a power flow diagram and show the direction of the line 
flows. 

(c) Check the power flow solution using the Ifgauss and other required pro- 
grams. (Refer to Example 6.9.) 


X t = j 0.0125 


Xt = j0.16 


jO.25 


© 

Op 

h - 4 

i 

1.25:1 2 


FIGURE 6.22 

One-line diagram for Problem 6.9. 
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6 . 9 . The one-line diagram of a four-bus power system is as shown in Figure 6.22. 
Reactances are given in per unit on a common MVA base. Transformers 7\ 
and T 2 have tap settings of 0.8:1, and 1.25:1 respectively. Obtain the bus 
admittance matrix. 

6 . 10 . In the two-bus system shown in Figure 6.23, bus 1 is a slack bus with V\ = 
1.0Z0° pu. A load of 150 MW and 50 Mvar is taken from bus 2. The line 
admittance is yi 2 = 10Z— 73.74° pu on a base of 100 MVA. The expression 
for real and reactive power at bus 2 is given by 

P 2 = 10| V^| | Vi | cos(106.26° -82 + 61 ) + 10|V 2 | 2 cos(— 73.74°) 

Q 2 = — 10jV^||Vi| sin(106.26° -S 2 + S{) - 10|V 2 | 2 sin(-73.74°) 

Using Newton-Raphson method, obtain the voltage magnitude and phase 
angle of bus 2. Start with an initial estimate of |V 2 |^ = 1.0 pu and S 2 ^ = 
0°. Perform two iterations. 



FIGURE 6.23 

One-line diagram for Problem 6.10. 


6 . 11 . In the two-bus system shown in Figure 6.24, bus 1 is a slack bus with V\ = 
1.0Z0° pu. A load of 100 MW and 50 Mvar is taken from bus 2. The line 
impedance is z \2 = 0.12 + y'0.16 pu on a base of 100 MVA. Using Newton- 
Raphson method, obtain the voltage magnitude and phase angle of bus 2. 
Start with an initial estimate of | V 2 1 ^ = 1.0 pu and = 0°. Perform 
two iterations. 



\ 

Z \2 = 0.12 + j'0.16 


1 V 

/ 

Vi = : 

.0Z0° 

1 > 


100 MW 
50 Mvar 


FIGURE 6.24 

One-line diagram for Problem 6.1 1. 
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6 . 12 . Figure 6.25 shows the one-line diagram of a simple three-bus power system 
with generation at buses 1 and 2. The voltage at bus 1 is V — l.OZO 0 per unit. 
Voltage magnitude at bus 2 is fixed at 1.05 pu with a real power generation 
of 400 MW. A load consisting of 500 MW and 400 Mvar is taken from bus 
3. Line admittances are marked in per unit on a 100 MVA base. For the 
purpose of hånd calculations, line resistances and line charging susceptances 
are neglected. 



MW Mvar 

FIGURE 6.25 

One-line diagram for Problem 6.12 

(a) Show that the expression for the real power at bus 2 and real and reactive 
power at bus 3 are 

P 2 = 40|V 2 ||Vi| cos(90° -S 2 + fc) + 20|V 2 ||V 3 | cos(90° - S 2 + <5 3 ) 

P 3 = 20|V 3 ||V 1 |cos(90° -Ss + h) + 20|V 3 ||F 2 | cos(90° - 5 3 + S 2 ) 

Q:i = -20(V 3 |Li| sm(90 o -<5 3 +5i)-20|V 3 ||F 2 |sin(90°-(5 3 +5 2 )+40|V' 3 | 2 

(b) Using Newton-Raphson method, start with the initial estimates of = 
1.0 -f jO and = 1.0 + jO, and keeping |V 2 | = 105 pu, determine the 
phasor values of V 2 and F 3 . Perform two iterations. 

(c) Check the power flow solution for Problem 6.12 using lfnewton and other 
required programs. Assume the regulated bus (bus # 2) reactive power limits 
are between 0 and 600 Mvar. 

6 . 13 . For Problem 6.12: 

(a) Obtain the power flow solution using the fast decoupled algorithm. Per- 
form two iterations. 

(b) Check the power flow solution for Problem 6.12 using decouple and 
other required programs. Assume the regulated bus (bus # 2) reactive power 
limits are between 0 and 600 Mvar. 

6 . 14 . The 26-bus power system network of an electric Utility company is shown 
in Figure 6.26 (page 256). Obtain the power flow solution by the following 
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methods: 

(a) Gauss-Seidel power flow (see Example 6.9). 

(b) Newton-Raphson power flow (see Example 6.1 1). 

(c) Fast decoupled power flow (see Example 6.13). 


The load data is as follows. 


LOAD DATA 

Bus 

Load 

Bus 

Load 

No. 

MW 

Mvar 

No. 

MW 

Mvar 

1 

51.0 

41.0 

14 

24.0 

12.0 

2 

22.0 

15.0 

15 

70.0 

31.0 

3 

64.0 

50.0 

16 

55.0 

27.0 

4 

25.0 

10.0 

17 

78.0 

38.0 

5 

50.0 

30.0 

18 

153.0 

67.0 

6 

76.0 

29.0 

19 

75.0 

15.0 

7 

0.0 

0.0 

20 

48.0 

27.0 

8 

0.0 

0.0 

21 

46.0 

23.0 

9 

89.0 

50.0 

22 

45.0 

22.0 

10 

0.0 

0.0 

23 

25.0 

12.0 

11 

25.0 

15.0 

24 

54.0 

27.0 

12 

89.0 

48.0 

25 

28.0 

13.0 

13 

31.0 

15.0 

26 

40.0 . 

20.0 


Voltage magnitude, generation Schedule, and the reactive power limits for 
the regulated buses are tabulated below. Bus 1 , whose voltage is specified as 
Vi = 1.025Z0 0 , is taken as the slack bus. 


GENERATION DATA 

Bus 

No. 

Voltage 

Mag. 

Generation 

MW 

Mvar Limits 
Min. Max. 

1 

1.025 




2 


79.0 

40.0 

250.0 

3 

1.025 

20.0 

40.0 

150.0 

4 


100.0 

40.0 

80.0 

5 

1.045 

300.0 

40.0 

160.0 

26 

1.015 

60.0 

15.0 

50.0 


The Mvar of the shunt capacitors installed at substations and the transformer 
tap settings are given below. 
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SHUNT CAPACITORS 

Bus No. 

Mvar 

1 

4.0 

4 

2.0 

5 

5.0 

6 

2.0 

11 

1.5 

12 

2.0 

15 

0.5 

19 

5.0 


TRANSFORMER TAP 

Designation 

TapSetting 

2- 3 


2-13 

mSå m 

3-13 

1.017 

4- 8 

1.050 

4-12 

1.050 

6-19 

0.950 

7- 9 

0.950 


The line and transformer data containing the series resistance and reactance 
in per unit and one-half the total capacitance in per unit susceptance on a 
100-MVA base are tabulated below. 


LINE AND TRANSFORMER DATA 


Bus 

No. 

Bus 

No. 

R, 

pu 

pu 

pu 

Bus 

No. 

Bus 

No. 

R, 

pu 

X, 

pu 

pu 

1 

2 

0.0005 

0.0048 

0.0300 

10 

22 

0.0069 

0.0298 

0.005 

1 

18 

0.0013 

0.0110 

0.0600 

11 

25 

0.0960 

0.2700 

0.010 

2 

3 

0.0014 

0.0513 

0.0500 

11 

26 

0.0165 

0.0970 

0.004 

2 

7 

0.0103 

0.0586 

0.0180 

12 

14 

0.0327 

0.0802 

0.000 

2 

8 

0.0074 

0.0321 

0.0390 

12 

15 

0.0180 

0.0598 

0.000 

2 

13 

0.0035 

0.0967 

0.0250 

13 

14 

0.0046 

0.0271 

0.001 

2 

26 

0.0323 

0.1967 

0.0000 

13 

15 

0.0116 

0.0610 

0.000 

3 

13 

0.0007 

0.0054 

0.0005 

13 

16 

0.0179 

0.0888 

0.001 

4 

8 

0.0008 

0.0240 

0.0001 

14 

15 

0.0069 

0.0382 

0.000 

4 

12 

0.0016 

0.0207 

0.0150 

15 

16 

0.0209 

0.0512 

0.000 

5 

6 

0.0069 

0.0300 

0.0990 

16 

17 

0.0990 

0.0600 

0.000 

6 

7 

0.0053 

0.0306 

0.0010 

16 

20 

0.0239 

0.0585 

0.000 

6 

11 

0.0097 

0.0570 

0.0001 

17 

18 

0.0032 

0.0600 

0.038 

6 

18 

0.0037 

0.0222 

0.0012 

17 

21 

0.2290 

0.4450 

0.000 

6 

19 

0.0035 

0.0660 

0.0450 

19 

23 

0.0300 

0.1310 

0.000 

6 

21 

0.0050 

0.0900 

0.0226 

19 

24 

0.0300 

0.1250 

0.002 

7 

8 

0.0012 

0.0069 

0.0001 

19 

25 

0.1190 

0.2249 

0.004 

7 

9 

0.0009 

0.0429 

0.0250 

20 

21 

0.0657 

0.1570 

0.000 

8 

12 

0.0020 

0.0180 

0.0200 

20 

22 

0.0150 

0.0366 

0.000 

9 

10 

0.0010 

0.0493 

0.0010 

21 

24 

0.0476 

0.1510 

0.000 

10 

12 

0.0024 

0.0132 

0.0100 

22 

23 

0.0290 

0.0990 

0.000 

10 

19 

0.0547 

0.2360 

0.0000 

22 

24 

0.0310 

0.0880 

0.000 

10 

20 

0.0066 

0.0160 

0.0010 

23 

25 

0.0987 

0.1168 

0.000 
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FIGURE 6.26 

One-line diagram for Problem 






